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A NONLINEAR TRAJECTORY COMMAND GENERATOR FOR A 


DIGITAL FLIGHT- CONTROL SYSTEM 
Luigi S. Cicolani and Stein Weissenberger* 
Ames Research Center 


SUMMARY 


A nonlinear Command Generator (CG) is synthesized for a digital auto- 
pilot capable of controlling VSTOL aircraft to track complex four-dimensional 
reference trajectories in an advanced terminal area air traffic control 
environment. The reference trajectory is assumed to consist of a sequence 
of straight-line or circular arc "legs" with discontinuities at the leg 
junction times. The CG generates smooth transitions from one leg to the next 
while maintaining a variety of operational constraints on maneuvering as 
dictated by safety, aircraft limitations, and passenger comfort, including 
bounds on velocities, accelerations, and acceleration rates. 

The basic approach regards the discontinuity as an error from the next 
leg of the reference trajectory which is then relaxed to zero in accordance 
with a suitable control law. A dual-mode control law is synthesized; it is 
globally asymptotically stable, has a maximal region of initial conditions 
using linear control, and uses minimum- time- to-origin control outside the 
linear region. Selectable control law parameters include the damping and 
settling time of the linear mode dynamics, saturation limits on velocity 
and acceleration excursions, and an acceleration rate limit. 

Transition maneuvers are generated independently along longitudinal, 
lateral, and normal flight-path axes, and the complete trajectory command is 
given by superposition of these transitions with the reference trajectory. 
Maneuver characteristics, such as duration, overshoot of the reference flight 
condition, and control activity, are optimized with respect to the time of 
initiating the transition from one leg to the next, and this optimization is 
done independently on each axis. 

Operational application of the CG was examined in detail in a simulation 
of a flight-control system with the augmentor wing jet STOL research aircraft; 
the basic repertoire of single-axis maneuvers and operational constraints 
are discussed, and the system behavior is tested on a rigorous STOL approach 
path and as affected by various approximations in the CG synthesis and types 
of disturbances found in the operational environment. The simulation results 
indicate that a satisfactory nonlinear system with general maneuvering 
capabilities was developed which satisfies the basic design objectives while 
maintaining a practicable degree of simplicity. 


^National Research Council Postdoctoral Research Associate. 


INTRODUCTION 


Problem Outline 

This report describes a trajectory Command Generator (CG) for an 
advanced digital flight-control system whose structure was develop,ed in 
reference 1. This structure provides a logical framework for designing auto- 
pilots for a wide range of aircraft, including CTOL, STOL, and VTOL, with 
capabilities of executing complex four-dimensional (one time and three space 
coordinates) terminal-area flight paths as directed by an air traffic control 
and guidance system or as selected by a pilot. Although the CG design is 
formulated here in the context of the structure of reference 1, its function 
could be embedded in most automatic flight-control systems independent of 
internal structure. Further, the structure of the command generator is 
independent of the aircraft being controlled. 

The control structure of reference 1 is outlined in figure 1. The 
Force Trimmap is the central element of the structure and computes the air- 
craft control command necessary to achieve an input commanded accelera- 

tion act . These computations are based on a model of the aerodynamic and 
engine forces of the aircraft being controlled (cf. ref. 2). The next level 
of the structure is the Trajectory Regulator (TR) , which produces compensatory 
acceleration commands on the basis of the error between the conmianded 

and estimated trajectories. The TR is intended to operate on relatively small 
path errors and hence is designed as a continuous- time linear, constant 
coefficient compensator, whose bandwidth is large compared to that of the 
command signals it receives from a higher level of the control structure. The 
combined Trimmap and TR form an acceleration controller structure that can 
execute and track commanded trajectories: 

|rc(t),vc(t),ac(t), < t < 

where r£ and V£ are three-dimensional inertial commanded position and 
velocity, provided these commands are kinematically consistent, within the 
aircraft flight envelope limits on velocity and acceleration, and provided 
that the trajectory errors remain sufficiently small that the total accelera- 
tion command act and actual velocity also satisfy these limits. 

It is the function of the trajectory Command Generator (CG) at the next 
level of the structure to ensure that these conditions are met. The CG is a 
feed-forward controller with position and velocity commands passed forward to 
the regulator and an acceleration command to the Trimmap. If trajectory 
errors become too large at, say, t^, the CG can initiate a transition maneuver 
that begins at the current estimated aircraft state ^(t^) , v(to) and relaxes 
the corresponding errors from the reference state to zero in accordance with 
a control law that bounds the velocity and acceleration excursions of the 
transition maneuver within specified limits. The synthesis of this control 
law is a central analytical problem of this paper. 
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The CG relieves the Trajectory Regulator and maintains acceptable 
commands in the following principal situations: (1) when very large tracking 

errors develop due to a period of large wind disturbances and it is desired 
to return to course; (2) when there has been a switch from one navigation 
sensor to another with a significant change in the estimated trajectory state 
and it is desired to reduce the newly sensed error from the reference; or 
(3) when a discontinuity in the input reference trajectory is imminent and it 
is desired to maneuver to the new reference trajectory. In this last situa- 
tion, the CG fills a role intermediate between the acceleration controller 
and the input reference trajectory. The source of the reference trajectory 
inputs may be a terminal area air traffic control and four-dimensional 
guidance system (ref. 3) or a trajectory tracking mode selection panel used 
by the pilot. These inputs are assumed to be a set of parameters sufficient 
to define the trajectory on the time interval [t^, t^] as a sequence of 
circular-arc or straight-line segments with constant path acceleration; each 
segment is defined from its initial position and velocity, the acceleration 
on the segment, and the corresponding initial time. A complete reference 
trajectory is then generated from a sequence S of such input parameters 


S 



M; t^ 





where v and refer to the speed rate and horizontal plane radius of 

curvature and are assumed constant on each segment [to^, input 

trajectory is assumed to be piecewise flyable; that is, each segment or leg 
corresponds to a steady or slowly varying flight condition within the flight 
envelope of the aircraft, but discontinuities in position, velocity, and 
acceleration can occur at each leg-switching time {tQ.}« The CG modifies 
the input trajectory to provide flyable transition maneuvering between 
successive legs by regarding the discontinuity as an error from the new leg 
of the input reference, which is then relaxed to zero in accordance with a 
suitable control law designed to satisfy operational constraints on aircraft 
maneuvering and provide command trajectories intermediate in bandwidth 
between the input commands and the trajectory regulator commands. The 
reference leg-switching times {to^} are usually known in advance so that a 
leg transition maneuver can also be initiated in advance of the reference 
switching time in order to minimize the required maneuver acceleration 
activity . 


Command Generator Structure and Transition Dynamics 

The structure of the CG is shown in figure 2. First, the position, 
velocity, and acceleration trajectory errors are formed in path axes: 


er = A (r - r*) 
p ps s s 

ev = A ^(v - V*) 

p ps^ s s 

ea^ = A^^(act - a*) 
p ps^ s s 


J 


( 1 ) 
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where (r*, v*, a*) and (r , v ) represent the reference and estimated 
trajectory states in a runway axis system with origin at the runway (fig. 3). 
The acceleration error is defined from act in equations (1) to obtain a 
desired continuity property in the output, as is noted below. These errors 
are transformed into the path axes system used by the transition dynamics; 
the path axes are sketched in figure 3 and are orthogonal with axes along the 
tangent of the commanded path (longitudinal axis) , normal to the tangent in 
the horizontal plane (lateral axis) , and normal to the tangent in the vertical 
plane (normal axis) . The transformation matrix from runway axes to path axes 
is given by 


A 

ps 


cos Y COS ^ COS Y sin m -sin y 
V V 


-sin ¥ 


cos ^ 


sin Y cos ^ sin ^ sin y cos y 

. V V 


( 2 ) 


where and y specify the heading and flight-path angle of the velocity 
vector (fig. 3). 


The error states (eqs. (1)) provide initial values for the transition 
dynamics states at those discrete times specified by the initialization 

switching logic 

e^Cto) = eVp(to) > (3) 

63 (to) = eap(t^) J 


The transition dynamics (TD) consists of three autonomous dynamic systems, 

l^(i) ^ j(i)(g(i)) ^ f(^)(0) = 0 , i = 1,2,3 (4) 


where the ith system governs the relaxation of position, velocity, and 
acceleration errors along the ith path axis, and the state vector for each 
system is 


e 


(i) 





(5) 


The complete maneuver command is constructed by combining these error 
relaxation maneuvers defined by the TD for each axis into the position, 
velocity, and acceleration vectors , e^ 2 > then transforming to 

runway axes and summing with the reference trajectory 
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rcg(t) = r*(t) + Ajge^(t) 
vcg(t) = v*(t) + 
acs(t) = a*(t) + A^^egCt) 


( 6 ) 


Coriolis accelerations associated with transitions can be neglected, as 
in equations (6), for the passenger operations considered here. These 
accelerations depend on the combination of path axis angular velocity (due 
to or y) transition excursions; however, rotation rates are held to 

small values in passenger operations, the effect is of short duration, and the 
size of the effect can be bounded by bounding the size of the initial errors 
that are to be relaxed by the TD. In summary, the CG structure in figure 2 
constructs a maneuver command that begins at t^ and at the then-current 
position and velocity, and acceleration command 


ic(t*) = ^(to) 
vc(t*) = i(to) 

= act (tp) 


> 


(7) 


and relaxes these states, by means of the TD, to the reference trajectory, 

\ / r*(t)> 

lim I vc(t) 

,ac(t) 



( 8 ) 


Path axes are used in the TD because aircraft maneuvering can be 
specified independently for each axis because of the nearly decoupled control 
that conventional aircraft have over motion along these axes. The path axis 
components of acceleration are related to kinematic variables of interest in 
aircraft trajectories by 



P 


= V cos 





(9) 


so that accelerations along each of the three axes control, respectively, 
speed, heading, and flight-path angle. Aircraft maneuvering can then be 
constructed as a superposition of single-axis motions; for example, the 
input reference trajectory is typically a sequence of legs that may be in 
static equilibrium » constant radius turns (constant R^) » or constant 

speed rate flight (constant v) or combinations of the latter two. The tran- 
sition dynamics modify this input trajectory by adding bounded time-varying 
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accelerations along all three path axes as required for good execution of 
the reference input path. For passenger operations, only modest accelerations 
are permitted; typical limits are 


|a<‘>| 
‘ P- ' 

< 

0.1 g 

•N 


‘ P ‘ 

< 

0.35 g 


> (10) 

1 P ' 

< 

0.125 g 

J 



Outline of Development 

The main analytical tasks of this work are the synthesis of the transi- 
tion dynamics (eqs. (4)) and the development of the TD initialization switch- 
ing logic. The transition dynamics are necessarily nonlinear because the 
trajectory is required to satisfy flight-operational bounds on velocity, 
acceleration, and acceleration rate. A dual-mode control is used to 
synthesize these dynamics, giving a linear system for sufficiently small 
initial errors and a nonlinear system for larger errors. In addition to 
satisfying trajectory bounds, the transition dynamics are required to be 
globally asymptotically stable; that is, errors must relax to zero starting 
from any initial value. Further, the dynamic system should be parameterized 
to permit adjustment of maneuver characteristics. Our approach to achieve 
these objectives uses an idea first reported by Kuznetzov (ref. 4), who 
defined the region of linear control using as a Lyapunov function the 
optimal cost function from the linear optimal control problem. The method is 
adapted to the present problem by a modification due to Willems (ref. 5), 
which permits a significant and straightforward enlargement of this region 
of linear control, A similar goal of linear region enlargement was achieved 
by Zachary et al. (ref. 6) where the approach, while direct, was analytically 
less simple and less easily applied to a variety of systems than the present 
method. Additional devices in the synthesis of the control law include 
Jordan canonical transformations and phase-plane analysis. A general approach 
to the design of dual-mode controllers is developed in detail in appendix A, 
modified to impose velocity bounds in appendix B, and applied to the present 
TD synthesis in the next section. 

The simplest initialization switching logic for the TD is to initialize 
at those times {to^} at which discontinuities occur in the reference input 
trajectory; this results in a unique transition for any given pair of 
successive legs. However, it is possible to optimize the commanded maneuvers 
regarding control activity and overshoot of the reference flight condition 
by adjusting the transition initialization times from the reference times. 

Such optimization is also a necessary property of an operationally acceptable 
flight system and corresponds to good "piloting technique" in manual flight 
operation when advance warning of a leg junction is available (e.g., 
ref. 7). Any leg of the reference can be extrapolated beyond its duration 
in the reference trajectory so that the initial conditions (eqs. 7) 
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for the transition to the next leg in the sequence are defined for any time 
earlier or later than the reference switching times. Consequently, the 
initialization times can be varied from the {to^} to improve properties of the 
ensuing transition. General switching criteria and closed-form solutions for 
the corresponding optimum switching times are derived for the dual-mode TD in- 
the fourth section of this paper. The results are given for arbitrary leg- 
junction discontinuities and are applied independently on each path axis; 
that is, the three autonomous dynamic systems of the TD are initialized 
independently and the results for each axis superposed as in equations (6) to 
yield an output command with the properties (eqs. (7) and (8)) at each 
initialization and in steady state. An alternative approach (ref. 8) consid- 
ers two successive legs extended indefinitely in time, with the first leg 
defined on (t < t^) and the second on (t > t^) and then formulates the 
transition as the solution of a linear optimal tracking problem on (-«> < t < «>) . 
While the approach in reference 8 allows systematic trade-off of reference 
path tracking excursions with control effort, it requires two sets of 
transition dynamics and, based as it is on linear theory, fails to enforce 
the various state and control constraints necessary for operational accepta- 
bility of the CG. 

The application tasks of the report are to derive a Command Generator 
algorithm for use in an application of the flight control structure of 
figure 1 to the Ames experimental powered-lift aircraft, the augmentor wing 
jet STOL research aircraft (AWJSRA) , and to evaluate the algorithm in 
simulation tests. A computer implementation is outlined and parameter 
values are selected to impose the kinematic constraints of the assumed 
passenger operational domain and the limits of AWJSRA capabilities. The out- 
put commands for the basic repertoire of single-axis maneuvers are examined 
to determine the effectiveness of the nonlinear TD in imposing the desired 
constraints and of the TD initialization switching logic in optimizing the 
maneuvers, and to determine the effect of approximations in the theory and 
design. Finally, the execution of a STOL approach trajectory is examined 
with a simulation of the flight control system and the AWJSRA. The results 
illustrate the general capability of the CG and of the flight control system 
to generate and execute satisfactory multiaxis maneuver commands. The test 
examines the effects of such features of the operational environment as 
off-reference initial conditions due to uncompleted earlier maneuvering, 
sensor switching events, and both random and steady disturbances. These tests 
include the principal disturbances to which the operational system will be 
subject. Future evaluation and development of the CG and system logic with 
increasingly realistic simulations and with pilot evaluation is planned. 

The structure and even many parameter values of the CG element of the 
flight control system are independent of the specific aircraft for which it 
is designed. Thus, details of the AWJSRA are suppressed in most of the 
discussion in this report and the design and test conclusions are largely 
applicable to a class of aircraft (CTOL and STOL) and to a class of opera- 
tions (passenger) . 
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TRANSITION DYNAMICS 


Figure 4 is a block diagram of the transition dynamics consisting of 
three decoupled dynamic systems that govern the error dynamics along each of 
the flight-path axes (these dynamics are initialized in accordance with 

(3))* For a single axis, the system consists of a control law, 

which is the basic acceleration command that defines the maneuver, 

and a filter, + 1)^\ which produces a smooth acceleration, e^. The 

resulting position (e^), velocity (e^) , and the acceleration (e^) commands 
from the three axes are combined with the reference trajectory to produce the 
output commands to the regulator and Trimmap elements of the autopilot. 

The control law, u^^^(e^^^), developed in this section, is linear for 
small errors but saturates at limit values for large errors and is also 
designed to limit velocity excursion, e|^) . Parameters of the linear mode of 
the control law are selected to obtain the desired damping and response time, 
while the saturation limits and switching boundaries of the nonlinear mode 
are designed to impose operational limits on the acceleration and velocity 
excursions which can be used for maneuvering in the neighborhood of the 
reference trajectory state. Because the control can be discontinuous, 

a filter is introduced to obtain continuous acceleration commands ^ 3 ^^ • 
Further, the filter parameter, r^i), is selected to limit maneuver accelera- 
tion rates, in accordance with the rate limits of aircraft capabilities 

to vary engine and aerodynamic forces. 

The third-order transition dynamics system for each axis seen in 
figure 4 will be used in the applications work. However, it is potentially 
useful to add a predictive modeling element that accounts for the predictable 
control lag effects on maneuver tracking; this expanded structure is outlined 
in the next section. The remainder of this section considers the control law 
design for a single axis, and the axis superscript ( ) is dropped from 
the notation. 


Transition Dynamics Design Objectives 


The various constraints to be imposed on the TD and the objectives of 
the control- law synthesis are defined next. First, the state equation for the 
TD is 


where 


e = Ae + bu(e) 


A = 


Ox 0 

0 0 X 

0 0 


= (0,0, t"^) 


( 11 ) 
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These dynamics are initialized to (1) reduce the regulator feedback error 
to zero, and (2) retain continuity of the total acceleration command to the 
Trimmap; that is, the components of _e are initialized as 

= ^p(to) - r*(to) 1 

^ 2 (^ 0 ) = Vp(to) - v*(tg) I (12) 

e^Cto) = [[actp(to)]] - a*(tQ> J 

It can be assumed that both the reference and current states in equations (12) 
are within the finite operational limits of the aircraft in velocity and 
acceleration. However, the acceleration command, act (t^) , might be outside 
the aircraft capabilities and should be bounded by the operational limits on 
acceleration (e.g., eqs. (10)) in order to ensure initialization at an 
admissible acceleration command. The bounded value is indicated by the 
bracket notation, [[ • ]] , which is defined as 

[[x]] = max{x ,mln{x,x }} (13) 

min max 

where ^ax understood to be known extreme values of x. The 

initial TD states are, therefore, contained in some bounded region denoted 
as 


63 ( 10 ) e [“® 30 m’® 30 M^ 


} 


(14) 


The transition dynamics generate the maneuver as an error from the 
reference state and then the CG command, referred to path axes, is 



(15) 


Constraints are imposed on the TD states based on operational considerations. 
These can be written in the form 


6^(1) e ^2 

63(1) e ?3 




43 d) 




i 


(16) 


These impose bounds on excursions about or overshoot of the reference 
trajectory; the velocity constraints limit the excursions in speed (or 
direction or flight-path angle) required to change position longitudinally 
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(or laterally or in altitude) and also limit any overshoot in maneuvering 
to a reference speed (or direction or flight-path angle). Similarly, the 
limits on e^ bound the acceleration overshoot in maneuvering to a new 
reference acceleration as, for example, when entering a turn. Lastly, 
bounds on the acceleration rate, e 3 , are imposed in order to limit the 
corresponding required aircraft control rates in accordance with the limits 
of aircraft control rate capabilities. 

The bounds ^ 2 ^ ^3 equations (16) are generally interior to the 

bounds on initial states in equations (14) , as illustrated in sketch (a) . 



Sketch (a).- Bounds on TD states. 

Thus, initial values of (e^, e^) can be outside 5^ and so that the TD 

should drive these states into ^2 5 3 in such cases. The bounds in 

equations (16) define the desired ranges of values in which the TD states are 
to be confined once entered; that is, the TD should have the property 


If 


^2 

then 

e2(t) e 52> t > to 

(17) 

If 


^3 

then 

® 3 (t) G ? 3 , -t > tQ J 


desired 

behavior 

for 

®3 

is obtained by bounding the control 

u as 


shown next. Solving the state equation for e^ yields 

e^Ct) = + ~ u(n)e'^'^^ drije (18) 

from which it follows that, if u is bounded by 

u(t) G (19) 


then ® 3 (t) satisfies 
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egCt) e[-e3n, + (esCto) + egje + (egCt^) - 63^)0 (20) 

Generally, if e^Ct^) is outside the bounds, then it is driven into and, 

once satisfied, the bounds remain satisfied for all further time. 

The rate constraints {^ 3 ^, also limit the acceleration rate of the 

CG trajectory in equation (15) , assuming that the reference command has 
constant path axis acceleration. Note that aircraft control rates required 
for the maneuvers are proportional to the commanded acceleration rates; hence 
the rate limits in equations (16) reflect aircraft limits on the rates at 
which engine and aerodynamic forces can be varied. The desired rate limits 
can be imposed by the selection of t . If u is bounded by equation (19), 
then it can be shown from equations ( 11 ) and (18) that the maneuver rate is 
bounded by 


T4g(t) £[-( 03 ^ + e3(t.,)),(e3^ - e3(t„))] 

Rate limits may then be imposed by requiring 

^ f®30m ®3m ®30M ®3M 

T > max < ; — , ; 

L ®3m ®3M 

Thus far, it has been shown that the structure of the transition 
dynamics permits maintenance of some of the constraints (eqs. (16)) provided 
that u and x are bounded according to equations (19) and (21) . It remains 
to indicate how velocity constraints are imposed and to complete the 
synthesis of the transition dynamics control law u(e^, e^, e^) . There are 
four important constraints and objectives that u must meet: ( 1 ) u 

( 2 ) the equilibrium of the resulting dynamic system (eqs. ( 11 )) must be 
asymptotically stable in the large, and (3) because of the simplicity and 
smoothness of linear control, the linear control region should be as large" 
as practicable. Furthermore, the time spent under nonlinear control should be 
minimized and there should be not more than a single transfer from nonlinear 
to linear control. Finally, (4) velocity excursion constraints must be 
imposed. 




( 21 ) 


Control Law Synthesis 


To accomplish the objectives listed above for the control law design, 
it is convenient to first transform the state space to Jordan canonical form 
(cf. ref. 9 ), changing from the system description of equations ( 11 ) to 


= ^2 1 

X2 = u(x^ .x^) J 


( 22 a) 
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and 



via the transformation equations 



(22b) 


(23) 


where the new state variables are 


(Xj, X2, X3). 


The inverse transformation is 


1x0 



0 1 


T 




0 0 1 . 


(24) 


Figure 5 is a block diagram of a TD structure employing these transforma- 
tions. Use of the canonical form greatly simplifies the design since the 
actual third-order system may now be controlled by a second-order law, 
u(x^, x^), as already indicated by the notation in equations (22a). This 
benefit is also retained when higher-order transition dynamics are treated. 
Although the control is now independent of x^, the characteristic root x 
for the Xo dynamics is selected by the designer and, more importantly, if 
the control is synthesized so that the (x^ , x^) subsystem is asymptotically 
stable, then the entire system is stable since u 0 implies ^ 0. 

A control law u(x^, x^) satisfying the four objectives listed above is 
derived in detail in appendices A and B and illustrated in figure 6(a). The 
state plane (x^ , x^) is divided into a region of linear control ^2^, which is 
compact and contains the origin, and the remainder of the state plane ^2^, 
on which the control is nonlinear. The linear control law is derived from 
quadratic optimal control theory and the region is computed from Lyapunov 

stability theory modified to obtain a large linear region with the property 
that all trajectories started in remain in it, that is. 

If x(1q) ^ then x(t) G t > t^ 

with 

u(x) ^ for all x G 

where 

X 5 (x^.x^) 
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The nonlinear control is designed for minimum time-to~origin from points in 
Qj^ within the control constraints. This differs from a minimum time-to-f^L 
control law but the distinction is not significant. 


The resulting control law and switching boundaries are summarized as 
follows: Let k = (k^, k^)^ be the gains of the linear control law. The 

region of linear control is given by the following combination of three 
regions of the state plane 


- {x:x'^Px < Vjj^} n [{x:-k'^x ^ ^2 ^ 

U {x:-k'^x ^ ^2 ^ 


(25) 


The switching parameter for the minimum time-to-origin saturated control is 
given as 


a = 


®3m^l + I ^2^ > ^2 ^ ° 

®3M^1 - 2 ^2^ ’ ^2 < 0 


(26) 


and the control law is 

/ 


u(x) = < 


-k*^x X ^ 

X G and [{a > 0,x, > -e,„} or {a = 0,x, > 0}] 


'3M 


2m 


X e and [{o > 0,x^ < -e^^} or {a < 0,x^ > 


2m 


X e and [{a < OjX^ ^ ^2M^ ^ ^ 


(27) 


The boundaries of £7^, the nonlinear region control switching lines (a = 0, 

^2 ^ "^2m> ^2M^ » control law u(x) are illustrated in figure 6(a). 

Some typical transitions that result from this control are shown in fig- 
ure 6(b). From an initial state in (at Xoi^ state is driven by the 

linear control to the origin without leaving Starting in state 

may be driven at saturated control into (from Xq 2 ) 5 into the switching 

line first (from Xq 3^* where the saturated control switches sign and the 
state is then driven into along the switching line; or into the velocity 

limit (from Xq^) where the control is simply nulled until the switching 
boundary is encountered; or, finally, the transient may begin in the region 
of saturated velocity and null control (at x^^), in which case the control 
remains null until the switching boundary is encountered. 


For the linear mode, the gains k of the control law, as well as the 
metric P and the parameter Vjjj used to define ^2^, can all be given in 
terms of the parameters of the characteristic roots of the maneuver dynamics, 
p, (damping ratio and time constant); suitable expressions, derived in 
appendix A, are 
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1 


p = 


^‘1 



CVJ 


\ - 24 + 6)max{e3^,e3^} f 


( 28 ) 


Thus, desired linear mode maneuver characteristics, such as critical 
damping, are easily imposed. For aircraft maneuvering, damping in the range, 
0.707 < ^ < 0.1, is desired while t^, the desired maneuver settling time, is 
based on passenger comfort restrictions or aircraft response limits on each 
axis. (The derivation of the linear control law is restricted by the condi- 
tion, p > 0.5, but this is not a significant restriction here.) In the useful 
special case of critical damping = 0.707) the control has some special 
properties. Assuming equal control bounds (e^^ = e 3 ]y[ = e^^) it can be sho\m 
that the end points of the pieces of the boundary enclosing have the 

simple coordinates shown in figure 6(c). From this, it is apparent that the 
linear control region expands with that is, if slower dynamics are 

specified, then larger initial state errors are relaxed with linear dynamics. 
In addition, the switching boundaries, a = 0, intersect those sides of the 
boundary of ^2^ on which the linear control saturates so that no discontinu- 
ity in u(t) occurs at the switch from nonlinear to linear modes for all 
those transitions which enter along the switching curve. 

The velocity excursion limits, imposed on the canonical 

velocity X 2 in this control law by the two regions of null control. Noting 
that x^ and u are equal by construction, the control law gives nonincreasing 
jx^l whenever is outside the limits; that is. 

If ^ $2 then ^ Ix^l < 0 (29) 

so that during a maneuver the state cannot violate these limits once they are 
satisfied; that is. 

If X 2 (to) e ^2 then X 2 <t) ^ ^ 


The maneuver may initialize at any value of x^ ; if X 2 (t^) is within the 
bounds it remains there, but if not, then either |x^| is reduced at the 
maximum possible rate or the state moves with constant velocity, x^(tQ), to 
the switching boundary in finite time, and then [x^i is reduced at the 
maximum rate and enters at some finite time, t^ 
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If X 2 (tQ) ^ ttisn x^Ct) ^ ^ 2 * 




(31) 


Thus, the canonical velocity is directly limited by the velocity bounds. 

The corresponding effect of the control law on the act\ial speed 
excursions, e^^ can be derived from equations (11) and (23), which yield 


= xe^ (32) 

from which it follows that e^ is always moving toward x^ 

If < x^, then 

If e^ > x^, then 

Then, as a result of equations (30) , (31) , and (33) , the actual speed also 
satisfies property (eq. (31)), except that the time at which e^ enters 
may differ from the time at which x^ enters. 

To summarize, the control law was synthesized to satisfy the four objec- 
tives listed at the start of the synthesis and the constraints (eqs. (16)). 

The synthesis utilized a Jordan canonical transformation (eq. (24)) which 
reduced the problem to one of stabilizing and controlling a second-order sub- 
system in (x^, X 2 ) (eqs. (22a)). In the resulting control law (eqs. (25)-(27)) 
the x-plane was divided into regions of linear and nonlinear control, 

In both regions the control u is bounded 


^2 > 0 


e2 < 0 


(33) 


— e ^ u 

3m 


e 


3M 


so that on any transition, c^(t) satisfies the same bounds in accordance with 
equation (20) ; the control is such that the velocity excursion constraints 
on X 2 (t) and C 2 (t) are satisfied in finite time from any point in the x-plane 
and, once satisfied they remain satisfied. The linear region is compact, 

contains the origin, is invariant (no trajectory can exit ^ 1 ,) 9 and the linear 
control provides asymptotic stability for the system from any point in as 

is shown in appendix A. Finally, for points not in the control provides 

minimum time to the origin and finite time to thus the equilibrium of the 

system (eqs. (11) or (22)) is globally asymptotically stable. 


Vertical Axis Transition Examples 

Results are now given that illustrate the control law and its application 
to vertical axis maneuvers to execute altitude and flight-path angle changes. 
The required maneuvering for these cases is almost entirely along the normal 
path axis; values of the transition dynamics parameters for this axis are 
taken as 
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®3m ®3M 


®2m ®2M 




T = 


0.125 g 
4.25 m/s 
(0.73,2) 
2 


Figure 7(a) shows the transition maneuver generated to execute a step drop 
of 60 m in the reference altitude for level flight; initial conditions are 

= (-60,0,0) 

The transition is shown as time histories of the output command (rc> ^ , 

vc|^^, ac|^)), the reference command (r*^^\ a*^^^), and the transition 

acceleration control u. In this example, the error states (e^^, e 2 , e 3 ) 
are the differences between the reference and output commands shown. As seen 
in figure 7(a), the aircraft descends to the new reference altitude without 
overshoot and at a descent rate that saturates at the speed limit, Thus, 

the velocity excursion limits for this axis act as the maximum descent and 
climb rates at which the CG executes altitude changes. The effectiveness of 
the velocity limiting logic in bounding both x^ and e^ is also seen in the 
state-plane plots of (x^ , x^) and (e^, e^) • The control history u(t) consists 
of a pushover period with saturation at ^ 3 ^, followed by a null period during 
which the descent rate is saturated at ending with a pullup period 

that is initially saturated at terminates in the linear control 

regime. Although u(t) is discontinuous, the corresponding filtered 
acceleration command ^^(t) and the output aCg(t) are continuous and rate- 
limited. The character of the entire maneuver is dominated by the nonlinear 
speed and acceleration saturation features of the control law, a phenomenon 
necessary to achieve acceptable maneuvering within operational limits even 
for the present modest altitude step example. A summary view of the CG 
function in this example is obtained by comparing the output command history 
(rcg, vcg, acg) with the input reference commands (r*, Vg, ag) in figure 7 (a). 

The same altitude step example with the addition of a large altitude 
rate tracking error of 6 m/sec and an acceleration error of 0.125 g (due to 
TR activity at the initial time) at the maneuver initialization time is shown 
in figure 7(b). The transition dynamics adapt to this error by initializing 
in accordance with equation ( 12 ) at 

e(to) = (-60,-6,0.125 g)^ 

This initialization nulls the apparent tracking errors seen by the trajectory 
regulator at tj and maintains continuity of the acceleration command to the 
Trimmap, act , within the bounds {e 3 Q^, output commands (vCg, aCg) , 

are now discontinuous at t^, and the early part of the maneuver differs in 
the behavior of (rCg, vCg) compared to figure 7(a) in order to relax the 
initial velocity error. 
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Figure 8(a) shows a flight-path angle change maneuver; the reference 
input path changes at to from level flight to descent at 4.3 m/sec. The 
initial conditions for the transition dynamics states for this case are 

e(t^) = (0,-4. 3, 0)T 

The commanded altitude history lags the reference path and eventually catches 
it without overshoot. The commanded descent rate necessarily overshoots the 
reference descent rate temporarily to attain the reference glide path. 

However, it is possible to eliminate the descent rate overshoot by initializing 
the transition to a new glide slope in advance of the reference switching 
tme as is commonly done in acquiring the ILS beam for final approach. 

Criteria and formulas for appropriately advancing the switching time are dis- 
cussed later in this report. The control u(t) in figure 8(a) corresponds to 
a pitchover saturated at e 3 ^ followed by a pullup to the reference flight- 
path angle saturated initially at This pullup and its associated 

reversals of the acceleration command and aircraft controls can also be 
eliminated by initiating the pitchover in advance of the reference switching 
time. 


Figure 8(b) compares the pitchover maneuver for five sets of transition 
dynamics parameter values listed in table 1 for the cases in figure 8(b): 


TABLE 1.- PARAMETER VALUES FOR EXAMPLE PITCHOVER MANEUVERS 



Case 

(P.Tl) 

X, 

sec 

m/ sec 

r ■ 

(1) 

Reference case 

0.73,2 

2 

4.3 

0.125 g 

(2) 

Slower linear dynamics 

.73,6 

2 

4.3 

.125 g 

(3) 

Reduced acceleration 
rate limit 

.73,2 

6 

4.3 

.125 g 

(4) 

Reduced speed limit 

.73,2 

2 

1.25 

.125 g 

(5) 

Reduced acceleration ' 

limit 

.72,2 

2 

4.3 

.06 g 


Reference case (1) is identical to that in figure 8(a). Case (2) shows the 
effect of increasing the time constant of the linear mode dynamics from 

2 to 6 sec. This expands the linear region 0,-^ considerably (as noted in 
fig. 6(c)), and the transient now occurs entirely in the linear mode. In the 
resulting maneuver, peak acceleration and velocity overshoot have decreased, 
but the maneuver duration has lengthened. Generally, increased x^ gives a 
smoother transient by expanding the linear mode region, but at the expense of 
increased maneuver time. 
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In case (3) , the time constant that limits acceleration rate t is 
increased to 6 sec. This expands the initial canonical state and the entire 
x-plane transition. Peak acceleration rate is markedly reduced as is peak 
acceleration, but a period of speed saturation appears, along with the largest 
speed overshoot among the cases shown; the position error state is also 
increased everywhere, and the maneuver is lengthened. In general, x is 
effective in controlling peak acceleration rates and should not be set higher 
than necessary to achieve this objective. 

In case (4), the speed excursion limit is reduced to 1.25 m/sec, and this 
is seen to be effective in reducing the velocity overshoot, but with a 
lengthened maneuver. 

Acceleration limits are reduced in case (5), effectively reducing the 
peak acceleration. The entire maneuver is slowed down and position error is 
increased everywhere due to the lower acceleration limit. These limits should 
thus be selected solely to control the peak acceleration. 


MANEUVER TRACKING PROPERTIES OF THE SYSTEM 


In this section the transition dynamics are considered in the context of 
the complete flight-control system. Maneuver tracking errors are analyzed 
and the use of predictive modeling of the aircraft response in the transition 
dynamics in order to prevent the use of feedback for maneuver execution is 
considered. The initialization rules and control law for the expanded transi- 
tion dynamics are also given. 

The maneuver commands are executed and tracked by the acceleration con- 
troller shown in figure 1. Tracking errors develop during all maneuvers due to 
control response lag, and these errors activate the Trajectory Regulator feed- 
back in the same manner as path errors due to external disturbances. Since 
the feedback is added to the maneuver acceleration command, sufficiently large 
tracking errors due to control lags can cause the total acceleration command 
to the Trimmap to exceed the acceleration limits that were carefully imposed 
on the trajectory command. We consider here how to modify the transition 
dynamics design so that, in the absence of external disturbances, the correc- 
tive feedback acceleration command during maneuvering is nulled or minimized. 
One approach, explored in the next section and applied in the AWJSRA system, is 
to reduce or minimize the acceleration activity (and, hence, control lags) 
required to transition between any successive pair of trajectory legs. Another 
approach, explored here, is based on predicting the tracking errors from a 
model of the aircraft. If such a model is included in the transition dynamics 
then the position and velocity commands transmitted to the Trajectory Regulator 
are the predicted aircraft response to the feed-forward maneuver acceleration 
commands and the feedback commands excited by maneuvering are eliminated to the 
accuracy of the model. Furthermore, the bounds on excursion velocity and 
acceleration and on acceleration rates ban now be imposed on the predicted 
aircraft states as well as the commanded states. 
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The block diagram of figure 9 represents the complete system of figure 1 
for a single path axis (the superscript (i) indicating the axis is dropped 
from the notation throughout this section) . Estimation errors and Coriolis 
accelerations are neglected in the discussion. The element G(s) denotes the 
transfer function between the commanded and actual accelerations. For conven- 
tional aircraft, G(s) principally models the engine, roll, and pitch dynamics, 
respectively, for the longitudinal, lateral, and normal path axes. The 
steady-state gain of G(s) is taken as unity. The acceleration disturbance 
includes the inevitable error in the Trimmap’s model of the aircraft 
forces. This error is assumed constant or slowly varying here and, together 
with G(s), represents the combined Trimmap and aircraft elements of figure 1 
for a single axis. An additional random disturbance noise, a^, represents 
the effects of wind turbulence. The transition d 3 mamics have been expanded to 
include an approximate model of the plant dynamics G(s) and an estimate of 
the Trimmap mismatch, a^. The dimensions of G and G are assumed equal, 
ignoring possible dimensional mismatch, and denoted as m. In the modified 
transition dynamics 03 is the commanded maneuver acceleration while y^ , e 2 , 
e^ are now the predicted acceleration, velocity, and position excursion states 
of the aircraft. The third-order transition dynamics of the previous section 
are a special case of figure 5 with 

m = 0 , G(s) = I , ^d ~ ^ 

Finally, the linear feedback law of the TR is seen in figure 9 and, for 
generality, it includes acceleration error feedback. 

An inspection of figure 9 suggests that the tracking errors z^, z^, 
can be eliminated in the absence of external disturbances provided all 
modeling is perfect (6 = G, a^ = a^) ; the transition dynamics are properly 
initialized to match corresponding aircraft states; the regulator gains, k^ , 
k^ , k^ are chosen to stabilize the acceleration controller; and the reference 
acceleration, a* , is constant. This is actually the case, as is shown next. 


Maneuver Tracking 

To determine tracking properties, a complete set of state equations for 
the system of figure 5 is necessary. For the TD these can be given as 


= ^2 


^ 


Yl = 72 


Vl = 


m 


?„--E « 


i^l 




= (u - eg)/! 


(34) 
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and for the aircraft. 


r = V 
P P 

= yi + a. + a 


Yl = Y2 


\ ( 35 ) 


ym-i = y- 


m 


ID 


= -r “lyi + a, (a* + 63 - + kgZg) 


Here, G(s) and G(s) have been given state variable representations in phase- 
canonical coordinates {y±} and {y^}, and the model parameters {a^^, 
i = 1, . . m}. Tracking errors for the complete system are now defined as: 


^ I 

Zt = r + e. - r 
1 p 1 p 

^2 = Vp + - Vp 

^3 = + yi - yi - 

2, = y2 - y. 


z . = y - y 

m +2 m -^m 


- a 


(36) 


Differentiating equations (36) and using equations (34) and (35) gives the 
error state equation: 


z, = z. 


z^ = z. 


z„ = z. + a^^ - a . 


z. = z. 


I (37) 


W = ■“2^ - • • • - + ^2^2 + (^3 + 1)Z3) 

m 


■ ^d> " E (“i - “i)yi 
i=l 


from which the following tracking properties are obtained. 
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Assume: 


1 . 

2 . 

3. 


then 


Perfect modeling: “ {o.^} , 

Constant disturbance and reference acceleration: 

= a* = 0 , ® 

k^, stabilize equations (37) (the roots of the matrix 

0 
0 


A = 


-Uiki 


i- - ot^(l + k^) - . . 

are all in the left half of the complex plane) , 


-a 


m 


1. lim ^(t) = 0 

t^ CO 


2- [z.(tQ) = 0] => U(t) =0, t > t^^] 


(38) 


Under the stated conditions, the errors tend to zero and the aircraft asymptot- 
ically approaches the CG trajectory independent of initial conditions and the 
transition maneuver. Further, by initializing the TD to null the errors, 

_z (to ) 5 ideal tracking with nulled errors at all times is obtained, independent 
of the trajectory. Under ideal tracking, the feed-forward commands would then 
excite no regulator activity and operational bounds imposed on the CG trajec- 
tory would be satisfied by the actual trajectory of the aircraft as well. The 
TD initialization rules for this ideal tracking are, using equations (36), 


e (t ) = V (t ) - v*(0 
2 o p o p o 

y (t ) = a (t ) - a*(t+) 
•^1^ 0^ p^ O'' p'^ o'^ 


y (t ) = y (t ) 
■’^m o •'m o 


(39) 
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This provides initialization rules for all TD states except the acceleration 
command e^ which can be initialized for continuity of the total acceleration 
command to the Trimmap to obtain continuous commands to the aircraft controls: 

eg(t^j) = [[actCt^)]] + - a*(tQ> (40) 

Here, the bracket notation [[•]} indicates that the operational bounds on air- 
craft accelerations are imposed on e^ and supersede, if necessary, the 
continuity condition. 

A case of practical interest occurs if no estimate of a^ is used 
(a^ = 0 ) and a^ has a nonzero constant value, but otherwise the assumptions 
of equations (38) hold. Then; 



Here, the trajectory regulator compensates for any constant disturbance or 
Trimmap model error at the cost of a steady position tracking error. However, 
the steady-state regulator output is now 


Aa 


P t 


-a 


d 


and is itself an estimate of the disturbance. 


Control Law Synthesis 


In the synthesis of the control law for the expanded transition dynamics 
(eqs. (34)) the objectives are the same as for the simpler third-order system: 
to achieve asymptotic stability with good transient properties and to satisfy 
the constraints (eqs. (16)), repeated here as: 


e 

e 

e 


2 

3 

3 


^ [-e 
e [-e 
G [-e 


2m’ 

3m’ 

3m’ 


"2M 


'3M 


'SM 


] 

] 

] 


Fortunately, the same synthesis procedure used in the previous section can be 
applied to the more complex system (eqs. (34)), as is demonstrated next. 


The constraints on acceleration excursion and rate, e^, ^ 3 , are again 
satisfied by bounding the control u as in equation (19) and selecting the 
filter time constant t in accordance with equation (21) . These results 
occur because e^ in equations (34) satisfies the same state equation as 
does e^ in equations ( 11 ) . 
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For the control law S 3 mthesis it is again convenient to decompose the 
state equations into Jordan canonical form. Transforming to canonical state 
variables, we obtain 


X2 

X3 


m +2 

9 

^m+3 



0 


0 

0 


.-1 


X. 


X. 


m +2 


X 


m+3 


+ 


0 

1 


r-1 


(41) 


where (A, b) are of the Jordan form (ref. 9). Since G(s) is assumed 
asymptotically stable the characteristic roots of A are all in the left 
half plane and the unforced dynamics of (X 3 , . - Xj^-(- 2 ) stable. The 

structure of equation (41) contains three independent subsystems, each driven 
by u(t), so that to obtain asymptotic stability of the entire state it 
suffices to synthesize u for stability of the (x^ , x^) subsystem 
((x^, x^, u) (0, 0, 0)). Thus, the control can again be taken as a function 
of only two states, u(x;^, ^ 2 ) , with precisely the same result previously given 
by equations (25)- (27). The corresponding canonical structure of the transi- 
tion dynamics is shown in figure 10 . 


The velocity constraints are imposed by the control law on the canonical 
velocity X 2 which again satisfies these constraints in accordance with 
equation (31). The commanded maneuver, as distinguished from the predicted 
maneuver, is given by integrating the acceleration command to the Trimmap; 
that is, by the states (e* , e^ in the auxilliary equations: 


e' = e' 




= e. 


= 


- (u 

T 


- eg) 


(42) 


The states {ei, 62 , yi) in equations (34) are the predicted aircraft response 
to the command 03 . Since the system (eqs. (42)) is identical to equations (11) 
it follows that the commanded velocity excursion e^ satisfies the velocity 
constraints exactly as did e 2 in the previous section in accordance with 
equations (30) and (31) . The behavior of the predicted velocity 62 relative 
to X 2 is governed by the following transfer function (obtained using expres- 
sions for X 2 ,B 2 from eq. (41) and the block diagram of fig. 9) 


= G(s) 

X^Cs) TS+1 


(43) 
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During saturation of the canonical^ velocity , X 2 (t) is constant at either “^2111 
or e^M then, recalling that G(s) has a steady state value of unity, 

equation (43) yields the desired steady-state saturation of the predicted 
velocity, that is. 

If x^(t) is constant then lim ^ 2 ^^^ ^ ^2 

t-> “ 


Finally, for reference, the canonical system was derived (using ref. 9) 
for two useful example cases: 

(a) a(s) - 


(b) G(s) 


1 

(t^S + 1)2 


where t, and t^, t ^ 0. The results are listed in table 2 (shown on 

the following page) . 


TRANSITION INITIALIZATION SWITCHING CRITERIA 


Thus far, the transition maneuvering has been constructed on the basis 
of a single leg of the reference input trajectory. The discontinuities at 
any leg switch are viewed as initial errors with respect to the new leg, which 
are to be relaxed by the TD, beginning at each of the reference leg switch 
times m = 1, 2, . . . , M) . However, maneuver characteristics can often 

be improved by initiating the transition in advance of the reference leg 
switch time, and criteria for selecting the switching time are indicated next. 
The practical value of such advance initialization is already well known in 
manual flight control for such operations as intercepting the final glide 
slope just before landing as well as for the more complex maneuvers of refer- 
ence 8. We seek a generalization of these manual flight tactics, suitable for 
our automatic flight-control system and applicable to the broader class of 
input commands considered here. 


It is useful to adopt notation which distinguishes the sequence of legs 
from which the reference trajectory is defined. This sequence is defined by a 
set of input parameters 


S 



.ro 


(m) _(m) .(m) 


.VO 



m = 1, 



(44) 


The elements of S are the states of leg-m at time • The trajectory for 

leg m can be defined at all times and is denoted (r^™^(t), v^™^(t), a(™^(t)). 

The reference trajectory switches from leg m to leg m + 1 at time and 

is thus defined on some flight duration [tO), tp] by 
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t 


E 



,t 


(m+1) 



m = 1, 


where = tr 



The leg-switching problem is illustrated in figure 11 for the pitchover 
example; the altitude histories for two successive legs of the reference tra- 
jectory are shown, with the reference switch time corresponding to the 

time of equal altitude for the two legs; the two legs are visualized as 
extending indefinitely far backward and forward in time. The relative error 
between the current trajectory and the next leg is defined at all times by 


|[^(t) - v(t) - v^^"^^^(t), a(t) - a^^"^^^(t)]; -oo < t < (46) 

Thus a set of error states is defined from which the transition can be 
initiated, and the initialization time can be selected from this set to obtain 
the most favorable initial state for the subsequent transient. The possibil- 
ity that the extension of leg-switching times to all t will bring the switch 
into conflict with other legs of the reference can be neglected here since 
in practice the parameters of S are selected such that the are close 

to the optimum switching times compared to the leg durations. Figure 11 shows 
three transitions generated by the TD previously defined, assuming perfect 
tracking and estimation at the time of initialization. Trajectory (a), which 
initializes at requires significant position excursion from the 

reference and, more significantly, a flight-path angle overshoot to rendezvous 
with the reference trajectory and maximum initial jerk. In contrast, 
trajectory (b) , initiated at time T prior to t^, reduces the tracking 
excursions, eliminates the velocity overshoot, and reduces the implied accel- 
eration excursions and acceleration rates for the maneuver, compared to (a). 
Finally, trajectory (c) is initiated too early and the TD attempts to climb to 
the backward extrapolation of the next leg. As in case (a), this requires 
maximum initial jerk followed by a control reversal and also the flight-path 
angle excursion exceeds the initial error that was to be relaxed. Thus, the 
initialization time generally has a significant effect on the characteristics 
of the ensuing transition synthesized by the CG. Criteria for the selection 
of this time are developed next. 


Initialization Time Criteria 

In this section, criteria are established to optimize leg-switch transi- 
tions with respect to the possible initial error states for the TD defined in 
equation (46). At any time t these states are defined by 
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(47a) 


Ei(t) = Ap3(t)(^3(t) - r(“+'>(t) 

A(t) = ApgCtXa^Ct) - a(”+^>(t) 

where Apg is the transformation to path axes from runway axes (eq. (2)). 
Then, for a single path axis, the corresponding initial values of the 
canonical TD states are, from equation (24) 


(i) 

(i) 


(t) 

(t) 


= E^^^(t) + x^^^E^^^(t) J 


i = 1, 2, 3 


(47b) 


The natural acceleration error A, in equations (47a) , governs the local 
evolution of the states E^, E 2 , and is distinguished from the constructed 
acceleration state, E^, used to initialize the TD. For simplicity, the formula 
for the canonical states above assumes the third-order transition 

dynamics although the criteria established here are valid for the higher 
order TD of the previous section. 

The criteria for selecting a best initialization time T for a single 
path axis depend on characteristics of the ensuing transition maneuver and 
differ for the nonlinear and linear modes. For the nonlinear mode transitions, 
the control u(x(t)) may undergo a discontinuous change with a sign reversal 
during the transition. The resulting commanded acceleration rate magnitude or 
jerk |e^^)| of the transition and the corresponding aircraft control rates 
(i.e., engine thrust, roll-angle, and pitch-angle command rates) are a local 
maximum at such a discontinuity. In addition, control sign reversals imply 
the possibility of velocity reversals and unnecessarily large velocity and 
acceleration excursions. Thus, the smoothness of the nonlinear mode transient 
strongly depends on whether there is a control sign reversal during the tran- 
sient and, consequently, a good criterion for initialization time selection is 
simply to prevent such reversals; that is, select T so that 

(j[x(t)] = 0 and x(t) ^ (48) 


where o is given by equation (26) and x(t) is given from equations (47) . 

This criterion ignores the character of the final linear mode portion of the 
transient since maneuver characteristics of interest are dominated by the 
nonlinear portion. Further, the locus of possible initial states x(t) defined 
in equations (47) may not intersect the switching curve, but this occurs only 
in cases such that the nominal transition is independent of the initialization 
time. 
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The effects of this choice of initialization time are illustrated in 
figure 12(a) for the pitchover maneuver. Maneuvers initialized at 0, 4, 
and 10 sec prior to the reference leg switch time are compared, assuming no 
errors in tracking the reference path at initialization. Reference time 
initialization (a) shows control sign reversal in u(t) and e^Ct) and velocity 
overshoot. Case (b) , which initializes 4 sec earlier and for which there is 
no control sign reversal, is an obvious improvement over (a); the velocity 
excursion for the maneuver does not exceed the discontinuity on the reference 
path, and the period at saturated control, maneuver duration, and peak accel- 
eration rates are all considerably reduced. Finally, (c) begins too early and 
introduces control sign and velocity reversals early in the transient where 
it attempts to climb to the backward extrapolation of the next leg. These 
three transitions are further compared in the x-plane in figure 12(b). The 
switching curves and boundaries separating regions of linear and nonlinear 
transitions are shown as dashed lines, and the locus of possible initial 
conditions for the transition, defined from equation (46), is given by the 
line L. Only nonlinear transitions can occur in this example because L never 
enters ^2^. The reference time initialization (at a) shows x^ overshooting 
zero until the switching curve, a = 0, is reached; a control sign reversal 
and reduction of the canonical speed follows. For the earliest initializa- 
tion (c) , the speed |x 2 | increases until the switching boundary is encountered, 
followed by control reversal and speed reduction. In both cases, the speed 
excursion on the transient, (x2)max “ (^ 2 ^niin» greatly exceeds the initial 
speed discontinuity to be relaxed. Finally, (b) , which initializes when L 
encounters the switching curve, enters the linear region without control sign 
or speed reversals and also requires the least maneuver time, time at 
saturated control, and speed excursion among these cases. Thus, the transition 
in this example is optimized (in the above sense) by initializing at the con- 
trol switching curve. 

For transitions initialized in the linear region, the initialization time 
can be selected to eliminate or minimize the magnitude of initial jerk to 
obtain the (locally) smoothest maneuver. Some consequences of this choice in 
the pitchover example are seen in figure 12(c). Here the locus L of possible 
initial conditions crosses the linear mode region. The zero-initial-jerk 
transition begins at (d) and uses the lowest velocity and position excursions 
and peak acceleration command compared to earlier initialization (at e) or 
later initialization (at a, c, or the reference leg switching time, b) and 
all nonlinear transitions are poorer in these measures than some linear transi- 
tions that can be initialized from L. Consequently, minimum initial jerk is 
the criterion used to initialize linear transitions; that is, the linear-mode 
jerk associated with initializing at t is given from equations (47) as 

|-x^(t)k*^^^ - (t) I 


and then T is selected such that 

J^(t) is minimized and x(T)^f^j^ (49) 
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An alternative linear-mode initialization criterion, studied in reference 8, 
would minimize the cost function of the linear optimal control with respect 
to the set of possible initial conditions. Criteria have been stated for the 
initialization time for both linear and nonlinear modes of the TD. By 
applying these criteria, transitions are smoothed, eliminating either initial 
acceleration rate or unnecessary control sign reversals. Such choices of 
initialization time also approximately minimize the extremes of acceleration 
and velocity excursions as well as, roughly, the overall acceleration activity 
and maneuver duration required, compared to other initialization times. The 
criteria given here lead to closed-form results for the initialization time 
which can be incorporated in the CG structure of figure 2 as the initializa- 
tion switching algorithm; these are derived next. 

The above criteria (eqs. (48) and (49)) apply independently on each path 
axis and, in general, distinct optimum switching times (T(i), i = 1, 2, 3) are 
obtained. The details of synthesizing the 3-axis command during the period 
between initializing the transition on one axis and on all three axes are 
given last in this section. 


Linear Mode Transition Initialization Time 

Throughout the derivation of the initialization time formulas, only a 
single axis is considered and the axis superscript is dropped from the 
notation. Third-order transition dynamics are assumed (G(s) = 1) and the 
formulas are derived for the time-to-go, AT, to the optimum initialization 
time for the convenience of on-line use of these results in the CG. If t is 
the current time, then optimum initialization time T is 

T = t + AT (50) 

For linear transitions, AT is calculated to minimize the jerk that 
occurs immediately following initialization; that is, to minimize over all 
At, the quantity 


l~k'^x(t + At) - E (t + At) I 

Jo(t + At) = ^ (51) 

T 

To evaluate the dependence of terms in equation (51) on At, let (E^ (t) , 
E 2 (t), E 3 (t), A(t)) be the current error from the next leg and let the future 
behavior of this error state be predicted assuming constant relative accel- 
eration, A(t) : 



+ 

+ 

+ 



1 

0 

0 


At At^/2 
1 At 
0 1 


( EjCt) 

E2(t) 

A(t) 


(52) 
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and constant relative acceleration coimnand, EsCt + At) = E 3 (t). This assump- 
tion is valid nominally (no disturbances, perfect tracking) and locally valid 
otherwise. 


The evolution of the corresponding canonical coordinates and initial 
TD acceleration and control are given from equation (52) using 


Xj (t + At) 
x^(t + At) 
eg(t + At) 
u(t + At) 


= E^(t + At) + TE2(t + At) 

= E^(t + At) + xE3(t) 

= E3(t) 

= -k^x^(t + At) - k2X^(t 4- At) 




J 


(53) 


The time-to-go, AT, which nulls initial jerk, is obtained by substituting 
equations (52) and (53) in (51) and solving 


0 = k3E^(t) + (k^ + Tk3)E2(t) + (xk^ + l)E3(t) 

[k^E^(t) + (k^ + xk^)A(t)]At + k^A(t)At^ 

+ - ■ ■ (54) 


In the remaining derivation the variable t is dropped from E^(t), E 2 (t), 
etc. when denoting conditions at the current time. The solution of equa- 
tion (54) yields the following cases: 


where 


AT 




undefined 

If 

A,E2 

-a 

If 

A = 0 

-b +Wmax{0,b^ - c} 

If 

A 0 



(55) 
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Three cases are distinguished in equation (55). In the first case 
(E^, A = 0) only a position step discontinuity occurs between the current 
state and the next leg as, for example, when a step change in altitude or 
lateral or longitudinal position occurs at the reference leg switch and there 
is perfect tracking. In this case, the initial conditions for the TD are 
fixed and independent of the initialization time; the initialization time 
is selected, therefore, on some basis other than minimum initial jerk. 

In the second case (E 2 ^ 0, A = 0) a velocity and, possibly, a position 
discontinuity, occurs but no acceleration change as, for example, where a 
discontinuity in speed or flight-path angle or direction occurs at the 
reference leg switching time and there is perfect tracking (e.g., fig. 12(c)). 
If, in addition, the reference legs are position-continuous at t^(Ej(tQ) = 0) , 
as is often the case, then equation (55) gives the optimum initialization time 
as 



This value is always earlier than ("leads") the reference time and is 
independent of E 2 (t^). This lead can be Incorporated in the reference 
trajectory itself by specifying a position discontinuity at t^ such that T 
and t^ are identical; from equation (55) the required discontinuity is 

Ei(t„) - 

In the last case in equation (55) an acceleration discontinuity occurs 
on the reference trajectory at to as, for example, when one of the legs is a 
turn or executes a constant acceleration speed change. Here, the quadratic 
equation (54) has two solutions which are either real or complex. They are 
complex if there is no time At at which the control, u(t + At) given by 
equations (53) can null the jerk. In this case, jerk is minimized by initial- 
izing at t = -b and. consequently this is taken as the optimum initialization 
time. 


If two real solutions occur, initial jerk is nulled at both times, but 
it can be shown that the acceleration magnitude |e 3 | increases following 
initialization at the earlier time and decreases following the later initiali- 
zation. Thus, the earlier solution produces a transition on which the 
acceleration excursion e^^^^ - ^3min necessarily exceeds the initial 
acceleration discontinuity E^ between the two legs; consequently the later 
initialization has been selected in equation (55) . 

Frequently, reference trajectories with an acceleration discontinutiy 
have continuous position and velocity at the reference leg switch time (E^ (t^) , 
E 2 (to) = 0) and then equation (55) gives: 
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Thus, the optimum time leads the reference leg switch time by an amount that 
is independent of the magnitude of E^Ct^). 

For on-line use of equation (55) it should be noted that if AT is 
negative, then the optimum initialization time has been passed and the initial 
jerk is an increasing function of time; consequently, jerk is minimized by 
initializing the transition immediately, and equation (50) is modified to 

T = t + max{0,AT} (58a) 

Further, the formulation in equation (55) is computationally sensitive at the 
boundaries between the three cases to the small nonzero values of or A 

which occur as a result of tracking errors and disturbances. If (A, E^) are 
zero on the reference trajectory, then small tracking errors have a large 
effect on AT but the resulting transition maneuver remains insensitive to 
the initialization time. If the reference trajectory is a velocity-step case 
then the optimum time is not sensitive to small nonzero values of A, but the 
formulation in equation (55) becomes ill-conditioned. Two devices are used 
in the applications work of this paper to avoid these two difficulties. 

First, the leg initialization time is limited to a short period, (t^ ” *^1 > 
t^ + "^2^’ neighborhood of the reference switch time so that 

equation (50) is further modified to 

T = max|^t^ " initi-|t + max{0,AT}, t^ + (58b) 

Secondly, the three cases in equation (55) are distinguished on the basis of 
approximate nulling of (A, E 2 ) rather than exact nulling by the use of thresh- 
old values, C2> ^3* conditions for each case are therefore: 


Case 

1 

If 

|A| 

< 

^3’ 

lEjl < ^2 'I 


Case 

2 

If 

|A| 

< 


IB 2 I > .2 

> (59) 

Case 

3 

If 

|A| 

> 

^3 

J 



Nonlinear Mode Transition Initialization Time 

For nonlinear transitions, initialization is taken at the intersection of 
the locus of possible initial states with the control switching boundary in 
accordance with the criterion given in equation (48) . This prevents control 
sign reversals at saturated levels. The intersection occurs at AT such that 
a[x(t + AT)] is nulled, where, as defined earlier, 


a(x) 


X 


X„ > 0 

3m 1 

2 2 

2 

, X 

1 2 

— -r- X„'^ 

c 

V 

>< 

3M 1 

2 2 

z 


( 60 ) 


Different functions define a(x) depending on whether the state x(t) is 
above or below the x^-axis and these are distinguished in the notation as 
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Equations (52) and (53) are substituted in equation (60) to 
evaluate the variation of o along the locus: 

0^‘*‘^(x(t + At)) = + I x^2) + [E^Ceg^ + A) + TA(e3^^ + E3)]At 

+ + A) ^ (61a) 

0^"^(x(t + At)) = (egj^x^ - I x^2^ + [E^Ce^j^ - A) + xACCgj^ - E3)]At 

+ A(e3j^ - A) (61b) 


The appropriate initialization time is selected from those times which null 
o in equations (61) and this selection results in the following cases: 


where 


r undefined 


AT 


• 

^-bi + -y^max{0,b^^ - c^} 


= 


u X, + 0.5 x„ X 
o 1 ' 2 ' 2 


u E^ 

O 2 


u 

O 2 


If A.E^ = 0 
If E 2 0, A = 0 
If A 0 


(62) 


-I 


'3m 


'3M 


x^ > 0 

x^ < 0 


b, = -7- 


o U + 0Eo 

a“ + 

A + |a| 


X. 


A - E 


+ T 


+ |Al 


c = 


2u^x^ + 0x^ 
A(u^ + A[) 


0 = sign (A) 




3m 

'3M 


A > 0 
A < 0 


Three cases occur in equation (62). First, if there is only a position 
discontinuity between the current state and the next leg (A, E^ = 0) the 
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initial conditions for the TD are fixed and independent of the initialization 
time, and no optimum time is defined. 


Secondly, if there is a velocity discontinuity but no acceleration change 
(E 2 0, A = 0) then the locus of initial conditions is parallel to the 

x^-axis and equations (61) yield a single time at which the locus intersects 
the switching boundary (fig. 12(b)). If the reference trajectory has continu- 
ous position at the reference leg-switching time and tracking errors are 
ignored, then the optimum initialization time always leads the reference 
switching “time as given by 


T = t - 
o 


T 


I ^2(^0) 

u 


(63) 


The required lead depends on the size of the velocity discontinuity |E 2 I’ 

A velocity change can be specified in the reference trajectory simply as a 
velocity step at some time t^, or, if very large lead time is required in 
equation (63), as a constant acceleration leg of duration E 2 /A beginning at 
t^. In either case the transition dynamics fill in the remaining details of 
the maneuver. The lead time in equation (63) can also be incorporated in the 
reference trajectory by specifying a position discontinuity at t^ such that 
T and t^ coincide, as given by 




(’ 4 ^) 


Finally, in the case of acceleration steps (A ^ 0) the locus of possible 
initial transition dynamics states is a parabola with axis of symmetry parallel 
to the x^-axis (see sketch (b)) and, from equations (52) and (53), can be 
given in the form 


^(X2(t + At) - x^^)2 
2 x, (t + At) - x 

1 IV 


with the time and location of the vertex given by 


At = 

V 




E 


3 


■ (>^1 - f (" * x ) > - ® 3 >) 

For positive (negative) values of A, the locus is traversed with increasing 
(decreasing) velocity; the x^-axis crossing occurs at time 


At 

X 
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Sketch (b).- Evolution of possible initial transition dynamics states; 

A > 0, 

and reaches a minimum (maximum) at the vertex. In perfect tracking, the 

vertex falls on the x.-axis. It is also noted that the switching curve 
(a(x) = 0 ) consists of two half-parabolas with vertexes at the origin and 
accelerations, and 03 ^^. 

The roots of equations (61) provide various possibilities for 0, 1, or 
2 points at which a(x) is nulled, as is illustrated in sketch (b) by the loci 
L^, L 2 , 1 * 3 , and their intersections with the two branches of the control 
switching curve. If A is positive, a single solution is selected in equa- 
tion (62) to minimize |a(+)|. Thus, if equation (61a) has two real roots 
(as for L 3 , in sketch (b)) the later one is selected and thus nulls a (“*"). 
If equation (61a) has complex roots (as for L 2 ) then the real part of 

the root is selected (At = -b^); this minimizes |cr and corresponds to a 
point on the locus between the x^- axis crossing and the vertex. 

These choices are based on the following rationale. If the locus inter- 
sects both branches (I< 3 ) , transitions initialized on either branch prevent 
reversals in u, but those initialized on the lower branch have increasing 
|e 3 | initially so that the acceleration excursion for the transition, 

^ 3 max ” ^ 3 min» necessarily exceeds the original discontinuity that was to be 
relaxed, E 3 . Transitions initialized on the upper branch decrease je^l 
at the maximum rate initially. If the locus intersects only the lower branch 
0 ^ 2)9 transitions initialized on the lower branch have slowly decreasing 1^31 
initially and the location and timing of the intersection is sensitive to 
tracking errors. In this and all cases for which no intersection of the upper 
branch occurs, initialization is taken at minimum This produces a 

reversal in u during transition, but the duration of transition, the time 
spent at saturated control, and the excursions in x^, x^ for the transition 
are all approximately minimized, and the initialization time is relatively 
insensitive to tracking errors. 

For negative values of A, an analogous rationale occurs by reversing 
the roles of the upper and lower branches of the switching curve, and the 
corresponding choice of roots of equation (61b) is contained in equation (62) . 
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If the reference trajectory has continuous position and velocity at 
the leg switch time (E^Ct^) = E 2 (tQ) = 0) and tracking errors are neglected, 
equation (62) gives 


■ Vui iVi 

Here, T always leads the reference leg switch time but by an amount less 
than T . 


(64) 


The same remarks noted for on-line use of the linear case results again 
apply to on-line use of equation (62) ; the transition is initialized as soon 
as the optimum time is passed (eq, (58a)), but the timing adjustment from 
the reference time is limited to a short period containing t^, as in 
equation (58b). Further, it is useful to distinguish the cases in equa- 
tion (62) on the basis of approximate nulling of A or E 2 as in equations (59). 


Synthesis of Three Axis Transition 

The preceding results define the initialization time for the transition 
along a single path axis. These apply straightforwardly in cases when the 
discontinuities on the reference trajectory are along only one path axis, as 
in the pitchover example. For maneuvers with discontinuities along two or 
three axes (e.g., simultaneous turn-entry, pitchover, and speed-change) to be 
relaxed by the transition, the initialization time optimization criteria 
apply to each axis independently and result in different optimum times for 
each axis. In this case, the TD can be initialized independently for each 
axis but the reference trajectory must also be switched to the next leg for 
that axis. The method of constructing this modified reference trajectory 
during the period between initializing one axis and all three axes is given 
next; the process consists of axis by axis superposition of the relative 
trajectory between the two successive legs. 

First, special formulas are given for extrapolating the next leg 
(leg m + 1) backward in time from the leg-junction at t . The discon- 

tinuity between the two legs at the reference switching time is given by 
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(m+l) _ . 
en; = A 
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^j.(m+l)J^^(m+l) ^j.(m+l)j 
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as referred to path axes of leg m calculated from v t ^ ^ - We define 
the difference between the two legs for t < to by assuming constant relative 
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acceleration in path axes en^, and using 
t(m+l); that is ' 


equations (65) as the difference at 


or. = or.^™'*'^^ 

en. 2 ^ ^ ' en ^ 


(m+l) 

^2 


en,(t) = enr'^' + enr'^^(t - 


(m+l), 

3 


( 66 ) 


enj(t) - en["+» + - t <“+»)+ | > (t - 


and the extrapolation of leg m+l backward in time from is 

constructed from the superposition of leg m and the relative trajectory 


.(m+l) 


(m+l) 


(m+l) 




(67) 


Equations (66) restrict the relative trajectory in path axes to one of 
constant relative acceleration during the period prior to t(^"^^) ; if both 
legs are straight lines with constant inertial acceleration, then the result 
in equation (67) is identical to extrapolating leg m+l from its input 
parameters* If either leg is a circular arc this is not the case, but 
the constant relative acceleration extrapolation is the appropriate one for 
use with the anticipation time formulas previously derived. 

The TD can be initialized at any time t < on the basis of errors 

from the extrapolated next leg given, using equations (47) and (67) , as 


^E^(t)\ 


E^(t) 


EjCt) 


.A<tV 





( 68 ) 


The first term in equation (68) is the difference between the estimated state 
and the state on leg m of the reference trajectory and is nonzero due either 
to tracking errors or to transitions which were initiated earlier but have not 
yet been completed. 
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Finally, if axis i is initialized at time the initial conditions 

for its transition dynamics are 

ej^^(0+) = j = 1,2,3 (69) 

and the reference path is switched to the next leg for axis i by means of 
the superposition 


where 


r*(t) = + Ajg(t)D(t)en^(t) 
v*(t) = Vg“^(t) + Apg(t)D(t)eii 2 (t) 
a*(t) = + Apg(t)D(t)eng(t) 


(70) 


D = 


<5i 



t < T 
t > T 


(i) 

(i) 


This reference trajectory coincides with the prior leg for times earlier 
than initialization of the first axis and with the next leg for times later 
than initialization of the third axis. 


A COMMAND GENERATOR ALGORITHM 


A Command Generator (CG) has been designed from the preceding theoretical 
analysis for use in an application of the control system structure of figure 1 
to Ames* experimental powered-lift aircraft, the augmentor wing jet STOL 
research aircraft (AWJSRA) . This design is defined next with flow charts and 
equations for the basic structure of the algorithm (fig. 13) and its principal 
elements; the reference trajectory generator (eqs. (73) and (75)), the transi- 
tion initialization algorithm (fig. 14), and the transition dynamics 
(fig. 15). A number of specializing choices occur in this application, and 
the rationale for these choices is noted in the discussion. The selection of 
parameter values is discussed in the next section. This algorithm has been 
subjected to and is the result of the simulation tests described in the 
following two sections. It is noted that the CG design structure, many 
parameter values, and the conclusions drawn from the simulation tests are 
independent of the specific aircraft being controlled within the class of CTOL 
and STOL aircraft, and for passenger operations; therefore, details of the 
AWJSRA are omitted in most of the discussion. 
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Input Parameters and Reference Trajectory Generator 


Figure 13 shows a flow chart of the CG logic executed each computation 
cycle in flight. The first element generates the reference trajectory for 
the current time and leg. The reference trajectory is assumed defined 
initially by a set of parameters which are the initial conditions for a set 
of M legs of the trajectory; 




, VO 


(m) 


• (m) 
r 


,R 


(m-l) 
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m = 1, 2, 


M} (71) 


These are the initial position and velocity in runway coordinates and the 
speed rate and horizontal plane radius of curvature for the leg. A number of 
properties are satisfied by these parameters and are assumed in the CG design. 
The set of initial leg times is ordered and strictly increasing and 

the reference path is intended to sequence through these legs switching from 
leg m - 1 to leg m at t^™) . The leg durations, normally 

longer than the transition settling time implied by the leg-junction discon- 
tinuities, (rol^"^^) - r^^^ (t , etc. , and the leg switching times are 
normally close to the optimum transition initialization times implied by the 
leg-junction discontinuities. These last two properties are not necessary 
conditions for the functioning of the CG logic of figure 13, but are expected 
properties of the class of operational inputs which have motivated the transi- 
tion dynamics concept and which underly the effectiveness of the transition 
dynamics initialization algorithm. It is also assumed that leg-junction 
discontinuities are limited in magnitude to avoid significant inconsistencies 
in the transition kinematics due to neglected Coriolis accelerations. Finally, 
the reference trajectory is assumed to satisfy operational constraints on the 
trajectory derived from passenger comfort considerations and the limits of 
aircraft capabilities. 


An analysis of the sources and definition of the reference trajectory is 
beyond our scope, but it can be generated by an ATC-4D guidance system, such 
as that described in reference 3. At the start of an approach to an airport, 
such a guidance system can define a reference path beginning at the current 
location of the aircraft and passing through a sequence of way-point positions 
(ro^^^) corresponding to the established approach routes of the airport; times 
and velocities at each way-point (to(^), vo|^)) and corresponding acceleration 
parameters (v(^), for each leg satisfy operational constraints for 

steady and nearly-steady flight conditions. 

The reference trajectory is generated from the input parameters listed 
in equation (71) . These parameters specify each leg as either a straight 
line with constant inertial acceleration or a helical path with constant 
radius of curvature. This set accommodates the restricted domain of passenger 
aircraft reference trajectories and of related guidance systems such as those 
of reference 3. For straight legs (Rq^ = 0), the acceleration vector is 
constant and aligned with the flight direction so that the reference trajec- 
tory is generated by the equations 
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where _S 3 ) are the runway axes. The horizontal plane coordinates are 

governed by a constant radius turn, where is positive (negative) for 

right-hand (left-hand) turns. It is convenient to define the unit tangent 
and radial vectors in the horizontal plane (u, m ) , the center of the circle, 
and the arc length from the initial position, s, in sketch (c) . The hori- 
zontal plane coordinates, (denoted x, y here) are given in terms of these 
quantities by 



The required quantities in equations (74) are calculated from the input 
parameters of equation (71) as follows: 
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Transition Initialization Algorithm 

The transition initialization algorithm in figure 13 is awakened and 
executed during a restricted period neighboring the reference leg-junction 
time, within this period the algorithm optimizes the time of initiat- 

ing the transition to the next leg and does so independently on each path 
axis. 


The timing adjustment period should be confined to the duration of the 
two legs being sequenced, but can be further confined in the present applica- 
tion to a short period, AT, on the order of 10 sec, in advance of t(^+^). 

This choice is based on the simplifying practice adopted here and elsewhere 
(refs. 3 and 8 ), of specifying reference trajectories with only single-state 
discontinuities at the leg junction (for each axis, only one of {E(i), e(i) , 

£ 3 ^^} is nonzero at The optimum initialization time for this class 

of reference trajectory discontinuities always precedes or "leads” 
by small amounts, as noted in equations (56), (57), (63), and (64), so that 
the restricted adjustment period accommodates this class and provides a margin 
for adaptation to off-reference conditions. The restriction on the choice of 
leg-junction discontinuities due to the above practice is not operationally 
significant and, in any case, the algorithm activates the transition at the 
best time in (t^^"^^) - AT, tC®"^^)) regardless of the discontinuities 
encountered. 

The transition initialization algorithm is shown in detail in figure 14 
and uses equations discussed in the preceding section. The algorithm, when 
awakened, predicts the optimum time to initiate the transition to leg m + 1 
on each axis, and whenever these times are reached the TD initialization indi- 
cator (ICTD) is turned on and will subsequently cause the TD to initialize. 

The internal indicator 6 (i) is also set for use in modifying the reference 
trajectory for this axis during the remainder of the adjustment period. If an 
axis is found to be ambiguous (|E 2 |j IE 3 I are below their threshold values so 
that the transition is negligibly dependent on its initialization time) , then 
another indicator, s(i), is used to postpone initialization until the algorithm 
turns itself off. Termination occurs if all three axes have been switched or 
are ambiguous (I(5(i) + s(i)) = 3) , or if the adjustment period has ended 
(t = ) , and consists of incrementing the leg counter m and turning off 

the internal indicators. A minor detail is the omission in equation ( 68 ) of 
the estimated acceleration tracking errors in evaluating the acceleration A, 
that governs the predicted locus of possible TD initial states. The omitted 
signal is dominated by random wind turbulence and is not usefully accommodated. 


Transition Dynamics 

The last element required for the CG is the transition d 3 mamics, which is 
shown in detail in figure 15 and uses equations previously discussed in the 
first two sections. Third-order transition dynamics with no modeling of plant 
dynamics are used; this choice was found adequate for the modest maneuver 
acceleration rates permitted in the passenger operations postulated for this 
application. The effects of this choice on system performance will be noted 
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in the simulation results. Initialization of the TD is controlled by the 
flag ICTD, which can be turned on either by the initialization algorithm 
described above or externally as a result of position sensor switching and 
associated large changes in sensor bias and in rg. The numerous input 
parameters of the TD 
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'3m 


3M’ m ’ 
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M . 


= 1,2,3 
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remain to be selected for each axis on the basis of operational considerations. 
Their selection is discussed in the next section along with illustrative 
transition examples for the basic single-state, single-axis leg-junction 
discontinuities. 

The saturation parameters • • •» denoted in figure 15 

simply as functions of the reference state to be calculated at the beginning 
of each leg, and appropriate equations are derived in the next section for the 
present design problem. This reflects the general nature of the operational 
constraints to be imposed in practice; they are limits on functions of the 
trajectory states from which corresponding saturation limits on the transition 
states can be calculated in terms of the reference flight condition. 


OPERATIONAL CONSTRAINTS, PARAMETER VALUES, AND BASIC TRANSITIONS 


The selection of parameter values for the AWSJRA application is discussed 
next in the light of kinematic details of the simplest and operationally most 
significant class of maneuvers; that is, transitions for single-axis, single- 
state leg- junction discontinuities. The TD saturation limits are used to 
enforce operational constraints on both steady state and transient values 
of various kinematic variables. These constraints may be generally applicable, 
such as those based on passenger comfort (ref. 8), or derived for the specific 
aircraft from safety margin considerations and aircraft performance capabili- 
ties (refs. 2 and 10). The list of constraints enforced below is considered 
preliminary. Considerable latitude for its expansion and elaboration is 
provided in the CG structure, however, in which the TD saturation limits can 
be calculated as functions of the reference flight condition and can be given 
different values for positive and negative excursions. 


Lateral Axis Parameters 

Lateral axis TD parameter values were selected as 
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For convenience, the notation ^ 3 e used to denote both minimum and 

maximum excursion limits when these are equal. In addition, parameters of 
the reference trajectory are indicated by the superscript ( )*. The velocity 
excursion limit, is given as a fraction of the reference aircraft speed 

at the initial time for the leg, v’^Ct^^^)), in order to limit the direction 
excursions, » used during transitions. These quantities are 

related by 

V* tanl'F^^ - ^*\ = |e^ 2 )| ^ ^( 2 ) 


and the selected fraction limits direction excursions below 14° at all speeds. 
During a sidestep, the aircraft departs no more than 14° from the reference 
flight direction and during a direction step the overshoot, if any, of the new 
direction, is limited to 14°. 

The lateral acceleration excursion limit ^ 3 ^^ is selected primarily 
to limit roll angle and heading rate commands appropriately. These quantities 
are related by 


^ * S.2 ^ ^ Y = (g cos Y “ ^ • £. 3 )^^^ (78) 

Desired operational limits on roll angle commands are 20° in steady state and 
30° transiently (ref. 8 ). Explicit heading rate limits are not precisely 
known, but 3.5°/sec is a typical maximum steady-state value for passenger 
operations. In the present context, equation (78) can be written in terms of 
reference and transition accelerations: 


(2) • ( 3 ) 

a* • £.2 + 63 = cos = (g cos y^ - ')tan 

.(3) 


(79) 


and Yq can be assumed a small angle while je^ | is limited to 0.125 g, as 
noted later in this section. Maximum roll angle and heading rate commands are 
then obtained as _ 

_l| ^an 


= tan 

' VC ' max 


1 

.(2) 


I 4' I ■ = I 1 + — ^ 

I VC 'max ' v' v_ 


( 80 ) 
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During sidestep and direction step transitions 

selected value of limits roll commands to 13° and Heading rates are 

below 3. 5° /sec for speeds above 65 knots. For acceleration steps it can be 
assumed that the reference trajectory satisfies steady-state limits, 

(taken as 20° and 3. 5° /sec here) or, equivalently, the reference acceleration 
satisfies: 


^ ^ ^ve ) 

Then, using the selected value of ^ 33 ^ limits transient roll angle commands 
below 33° in the worst case and maximum transient heading rate commands are 
below 6 . 6 ° /sec for speeds above 65 knots. Ordinarily, however, the TD initial 
time optimization effectively eliminates or minimizes overshoot of the 
steady-state roll angle and heading rate during acceleration steps and the 
margin for transient excursions is used by the TD only to adapt to off- 
reference initial conditions. 

The parameters > limit the initial absolute acceleration com- 

mand. Its selected value limits the initial roll command to 30°, approxi- 
mately, in accordance with equation (78) and the desired transient roll com- 
mand limits. 

The filter time constant is selected primarily to limit roll rate 

commands within the desired operational limit for passenger operations, 

10°/sec (ref. 8 ). Transitional jerk is given from equations (11) by 




|u<^> - a<^> 


,< 2 ) 


and, from the derivative of equation (79), the extreme roll rate command is 
given by 
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For position and direction step transitions, the maximum jerk and worst-case 
roll rate command occur at saturated control reversals: 


or 



2e 
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= 0.13 g/sec 



max 


9°/sec 


For turn entry or exit with good tracking, the maximum possible jerk and roll 
rate command occur initially in the transition 
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or 


< 


(g tan + € 3 ^^) = 0.19 g/sec 



T 




VC * max 


12 °/sec 


Finally, in the worst-case initial conditions (maximum turn from saturated 
opposite roll and simultaneous unfavorable normal axis maneuver commands) 
the jerk and roll rate are 

(g tan 30° + g tan = 0.38 g/sec 

T ^ 


N>. 


VC ' max 


22Vsec 


These estimates are, however, pessimistic because the TD initialization time 
optimization results in initial conditions for which the initial control, e^, 
reduces the initial jerk, and initial rates are below 5°/sec for turn entry 
and exit, and below 16° /sec in the worst case. 

The time constant scales the size of the linear mode region of 

the state plane as discussed previously. In the present work, is 

selected on each axis to be slower than the characteristic times of the plant 
dynamics in the simulation. This linear region is shown in figure 16 for the 
lateral axis. Sidestep transitions are initiated on the x^-axis; sidestep 
distances greater than 16 m require nonlinear transitions, such as those 
initiated at or x^^ and usually use maximum excursion roll angle 

commands and can encounter saturated control reversals (from Xq;^) with 
corresponding roll rate extremes under 10° /sec. Direction steps in excess 
of tan“^ /v*] require nonlinear transitions; the state follows loci, 

such as or L 2 in figure 16, until the switching boundary is encountered 

at which point the transition is initiated in accordance with the optimum 
initialization time algorithm and usually uses maximum excursion roll angle 
to turn to the new direction. Optimum initialization of direction steps up 
to tan~^[(AT - T)e^|)/v*] is accommodated by the advance period, AT, as 
determined from (63) . Finally, for the given parameter values, all accelera- 
tion step transitions are initiated in the linear region; the state follows 
a locus, such as L 3 in figure 16, until it enters after which the 

transition is initiated at the point of minimum initial jerk; this occurs at 
the boundary of in the common case that the acceleration step exceeds 

e(^). These acceleration step transitions require very little roll command 
overshoot and roll rates are well below 5° /sec as a result of the initiali- 
zation optimization. The lateral axis single-state reference trajectory 
discontinuities — position, velocity and acceleration steps — correspond, 
respectively, in flight operations to sidesteps to a parallel path, step 
changes in direction of straight flight paths, and entry to or exit from 
steady turns or reversals of turns. Lateral axis transitions result in 
rotations, of the path axes and related longitudinal axis cross-coupling 

activity in some of the basic transitions to be discussed next. 
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Sidestep Transitions 


Some kinematic variables from a typical sidestep to a parallel path at 
a distance of 250 m are shown in figure 17(a). As seen in the x-y trajectory 
plot, the TD generates a turn away from the reference direction, followed by 
a straight path to the parallel path and a return to the reference direction. 
Time histories of ^vc» ^vc» ^vc demonstrate the effectiveness of the TD in 
limiting these variables during the transition in accordance with the selected 
bounds. 

The sidestep distance for which the present TD can generate acceptable 
transitions is limited by several effects of neglected longitudinal axis 
cross coupling and Coriolis accelerations. First, a reduction in speed is 
commanded as a result of the combined lateral velocity excursion and path 
axis rotation, during the maneuver (see vc in fig. 17(a)). However, it 

is necessary to increase speed during the sidestep in order to maintain 
position along the reference trajectory, but no longitudinal axis transition 
is generated to provide the proper speed variations since the initial condi- 
tions for the longitudinal axis transition are zero. These effects must then 
be compensated by the control system feedback; that is, position errors and 
corresponding corrective feedback acceleration commands (A^ • p^) develop 
in order to increase speed (v) as required. This effect is proportional 
to the direction change and is limited by the TD bound on e^. 

Secondly, the TD excursion trajectory (Arcg, Avcg, Aacg) is kinematically 
incompatible longitudinally since the neglected Coriolis terms are nonzero 
and principally in the longitudinal direction in this case. A significant 
longitudinal position discrepancy develops at the rate which is 

greatest during the initial turn from the reference direction and dominates 
the feedback signal Aa. • p^ during this period; the discrepancy is 
proportional to sidestep distance so that the aircraft longitudinal accelera- 
tion capability is saturated by the feedback at sufficiently large sidestep 
distances. However, satisfactory sidesteps to very large distances can be 
obtained alternatively by incorporating in the reference trajectory the 
straight-line segment joining the two parallel paths; this reduces the transi- 
tions generated by the TD to direction steps with much reduced Coriolis 
effects. In any case, the use of feedback to accomplish trajectory commands 
is undesirable. 


Direction-Step Transitions 

The TD generates a turn to the new reference direction in this case, with 
initialization time selected according to equations (56) or (63) . Kinematic 
variables of interest from a typical direction change of 30'^ during level 
flight at 100 knots are shown in figure 17(b). Three cases are shown; they 
correspond to initialization at (1) the reference leg switch time, (2) the 
optimum initialization time, and (3) an intermediate time. The effective- 
ness of the initialization time optimization in reducing or eliminating 
excursions from the reference path, overshoot of the reference direction, roll 
command activity and reversals, roll rate peaks and activity, and maneuver 
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duration are apparent as well as the effectiveness of the TD in bounding 
peak values of roll rate, roll angle, and direction overshoot in accordance 
with the selected parameter values and for all cases shown. The algorithm 
allows initialization up to 10 sec in advance of the reference leg switch time, 
and this period allows for optimum initialization of direction steps up to 
tan”^ (20/v’^) . The 30°-step example uses most of this advance period. Larger 
direction changes can be accommodated with larger advance periods, but these 
transitions are better accomplished in the present application by including a 
constant radius turn in the reference input trajectory. 

A longitudinal axis transition is also generated for the direction step 
and has initial conditions corresponding to the apparent longitudinal velocity 
step at the leg-junction: 

= ( 0 ,v*(l - cos A¥^), 0 ) 

The discontinuity size increases with |a^vI significant size at 

30°. This transition is shown in figure 17(c) for the present example and 
significant acceleration and velocity activity is commanded. Since the trajec- 
tory velocity command results from the combined lateral and longitudinal 
transitions , 


VC = V* + A ^e^ 
s s ps 2 

the speed command variation, vc, (shown in fig. 17(c)) depends on the 
transition velocity coordinates of both axes. In principle, a direction step 
can be accomplished with no speed variation or corresponding longitudinal 
axis activity. In the present design the speed variation results from 
neglected Coriolis terms, whose effects are seen in the inconsistency of the 
variations of the actual and commanded speeds. 


Lateral Acceleration-Step Transitions 

Lateral axis acceleration steps occur at any change in the turn radius of 
the reference flight path, such as at entry to a steady turn from a straight 
flight path or at exit from a turn. Figure 18 shows results from entry to a 
level turn at 100 knots, with radius such that the maximum bank angle allowed 
for trajectory commands (20°) is required. The TD initial conditions corre- 
spond to the leg-junction trajectory discontinuity 

^ (0,0,0.36 g) 

The TD generates a brief transition which essentially smoothes out the roll 
step command of the reference path in accordance with the roll rate limits 
and adjusts the transition initialization time to avoid roll or direction 
overshoot in tracking the reference path (see fig. 18(a)). In this example, 
the transition acceleration command ^ 3 ^^ initializes outside the bounds 
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on control (±0.2 g) , but is relaxed to zero monotonically with very little 
overshoot of zero. 

Further details of this transition are shown in the x-plane plot 
(fig. 18(b)). The locus, L, of possible initial conditions for the turn 
entry, and transitions are shown, beginning from three points corresponding 
to (1) the reference leg junction time, t^, (2) the optimum time (t^ - 2), 
and (3) an earlier time (to - 6). The optimum transition begins where L 
enters and is entirely in the linear mode; the initial jerk is minimized 

at this point but is necessarily nonzero because is outside the bounds 

e(|) . The locus L has the property in this example that neither initial- 
ization criteria ((48), (49)) can be satisfied; that is, if t*, t^ are the 
optimum linear and nonlinear mode initialization times then 


x(t*) f and x(t^) $ 

For such cases the on-line initialization algorithm selects the earliest time 
that one of the following conditions is satisfied 

1. t > t* and x(t) e 

2. t > t* and x(t) ^ Q. 

NL L 

3. t > t 

o 

In this instance, (1) is satisfied at entry to and a minimum initial-jerk 

transition is generated with very little roll angle overshoot. The advance 
time required for optimum linear mode initialization of acceleration steps is 
independent of the size of the steps (see eq. (57)) and, for the lateral 
axis, is within the 10-sec advance period used in the algorithm. 

The effects of initialization time on the control (u) , on roll angle 
($Vc» position tracking excursion (e(^)) are also shown in 

figure 18(b). The transition initiated at the reference time is nonlinear; 
maximum negative control is used in an effort to "catch up" with the step roll 
change of the reference trajectory and corresponding roll rate command exceeds 
the desired 10° /sec limit. The roll command also overshoots the reference 
roll angle to the excursion limit and into the margin reserved for transient 
activity. 

In contrast, the optimum time initialization reduces roll rates well 
below the desired limits and virtually eliminates the roll overshoot; and, in 
comparison with other transitions in figure 18(b), it requires the least 
duration, control activity, roll overshoot and activity, and path excursions 
from the reference. 


Finally, it is noted that no longitudinal axis activity is generated 
during these acceleration steps either from the TD or from feedback due to 
neglected cross coupling. 
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Longitudinal Axis Transitions 


For the longitudinal axis, the single state reference trajectory 
discontinuities correspond respectively to position steps forward or backward 
along the flight path, step changes in airspeed, and speed rate steps. Their 
transitions are qualitatively the same as the corresponding position, velocity, 
and acceleration step transitions for the other axes and need be reviewed only 
briefly. Path axis rotations and related cross-coupling effects are absent. 

The following parameter values were selected for the AWJSRA: 


2m 

3m 
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-g sin 
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2 sec 


= 0.065 v*(to) 

= 0.06 g 
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(83) 


(p 


(l) 



(0.707,6 sec) 


These values reflect the limits of aircraft capabilities rather than 
passenger comfort restrictions. The velocity bound limits speed excursions 
used, for example, during position adjustments. This bound was selected as a 
fraction of the reference speed; it satisfies the small speed reduction margin 
available at minimum operating speeds and, at higher speeds, it bounds the 
implied aircraft control activity. 

The acceleration limits ^ 3 ^^ > selected in relation 

to the aircraft capability to generate longitudinal forces. These are the 
static engine and aerodynamic forces and they balance aircraft weight, as 
well as supply the trajectory acceleration, that is. 



where fc is the required applied acceleration for the trajectory and is 
to be supplied by the aircraft engine and aerodynamic forces. The longitudi- 
nal component, f^»p^, is denoted A^ and is given (in g’s) by 

sin Y + vc 

= _£ (84) 

c g 

Let (Au^, Au^} denote the aircraft limits on longitudinal specific force 
capability. These are determined from a variety of operational constraints on 
aircraft control usage and margin requirements and, in general, vary with the 
reference flight conditions (ref. 2). For passenger aircraft generally these 
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limits are modest and for the AWJSRA they are about ±0.15 g or less. It can 
be assumed that the reference flight conditions (y*> v*) provide a margin for 
transient longitudinal forces, on the order of 0.06 g, and this value is 
selected as the longitudinal acceleration excursion limit ^ 30 ^* However, it 
should be noted from equation (84) that this margin is shared by both normal 
and longitudinal axis variables, Yj v. Ordinarily, in velocity or acceleration 
steps, the initialization time is optimized to prevent overshoot of the refer- 
ence values of either v*, and the transient margin is used only in some 

cases of off-reference initial conditions. The initial longitudinal accelera- 
tion command is limited by equations (83) to reflect the 

approximate limits on the longitudinal specific force. 

The rate-limiting time constant x(^) is well above the time constant for 
the aircraft response to longitudinal acceleration commands, which are var- 
iously provided in the case of the AWJSRA by the engine, nozzle, and angle of 
attack variations depending on configuration. 

The region of linear control and switching boundaries for the longitudinal 
axis TD is shown in figure 19(a). Position steps in excess of 11 m require 
nonlinear transitions and maximum accelerations. The velocity bound is usually 
encountered; this bound increases with speed and intersects speeds 

below 109 knots. Speed steps in excess of 7 knots also require nonlinear 
transitions. 

Some results from a backward position step of 80 m during level flight at 
100 knots are shown in figures 19(a) and (b) . The position error ep) is 
relaxed from its initial value to zero monotonically by the TD with very little 
overshoot; the speed command vc is reduced to the value permitted by the 
excursion limit, where it is held until the position error is reduced suffi- 
ciently, and then returns to the reference speed with very little overshoot. 

The acceleration commands illustrate the effectiveness of the acceleration 
and acceleration rate limits of the nonlinear transition dynamics. 

A 20-knot speed-step transition during level flight is shown in fig- 
ures 19(a) and (c) . These results indicate the effectiveness of the TD 
initialization algorithm in eliminating saturated control reversal and corre- 
sponding overshoot of the new reference speed. It can be calculated that the 
10 -sec advance period provided in the present design allows optimum initial- 
ization of speed steps up to 19 knots. Larger speed changes are better 
accomplished by including a constant acceleration segment in the reference 
trajectory. In addition, the 10-sec advance period also accommodates all 
acceleration steps. 


Normal Axis 

Transitions corresponding to altitude and flight-path angle steps for 
the normal axis were discussed in a previous section. Transition dynamics 
parameter values selected for the AWJSRA are 
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The velocity bound is designed to limit flight path angle excursion amplitudes 
relative to the reference value to 3° and to limit altitude rate excursions 
to 6 m/sec. These are appropriate passenger operation limits on flight path 
angle and altitude rate during altitude step transitions. Further, it is 
assumed that values of y*, v^ on the reference path are restricted to pro- 
vide a margin for such flight-path angle excursions within the limits of 
aircraft longitudinal specific force capabilities, in accordance with 
equation (84) . 


The acceleration limit is a passenger comfort constraint (ref. 8 ) and is 
well below the minimum aircraft lift margin capability of 0.4 g to 0.7 g 
required for safety at all reference flight conditions. 


The region of linear control is the smallest for any of the three 
axes and is such that nonlinear transitions with saturated control are used 
for altitude steps above 2.5 m and altitude rate steps above 2.5 m/sec. 


The 10-sec advance period allows optimum transition initialization of 
all flight-path angle steps within the AWJSRA capabilities so that it is 
unnecessary to consider the use of any constant normal acceleration leg in 
the reference input trajectory. 


Finally, it is noted that a discussion of the important landing flare 
maneuver is omitted from this paper. The reference trajectory command is 
specially tailored in this case so that the transient trajectory response 
of the aircraft satisfies more detailed requirements than for the transitions 
so far discussed and more than are necessary elsewhere on terminal area flight 
paths . 


Remarks 

A task of this section has been to define the desired operational con- 
straints on trajectory commands and relate them to the excursion limits of the 
TD. For this purpose, various kinematic variables of concern (v, A^, $y, 

$y, Y, a • P 3 ) and constraints on their steady-state and transient values were 
introduced.” These are constraints on functions of the CG command and can be 
generalization to a list of the form 
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{gj^(rc,vc,ac,ac) > 0 , i = 1,2,. . . , N} 
which reflect operational practices. 

The input reference trajectory is assumed to satisfy the constraints on 
steady- state kinematics 

{gj^(^*,v*,a*,^) = 0 , i = 1,2,. . N} 

The constraints can also be given in terms of the reference trajectory and 
transition states 

{gi(r*,v*,a*; e^,e^,e^,e^) > 0 , i = 1,. . N} 

from which boundaries on the excursion states and corresponding saturation 
parameter values can be defined as functions of the reference flight condition. 

The TD parameters were demonstrated to be effective in general in 
imposing the desired linear mode transient characteristics and the desired 
bounds on commanded velocity and acceleration excursions about the reference 
flight condition. Further, the transition initialization time algorithm was 
demonstrated to be effective in approximately minimizing control activity 
and overshoot of the reference acceleration, velocity, and position during 
the transition- These characteristics are considered necessary conditions 
for an operationally acceptable trajectory command generator; they are 
features of good piloting technique which the present design attempts to 
generalize for use by an automatic flight-control system coupled to an Air 
Traffic Control/4-D Guidance system having a broad reference trajectory 
command set at its disposal. 

The use of the CG as part of the flight-control system of figure 1 and 
the system behavior during the execution of an approach trajectory and in the 
presence of off-nominal events and disturbances is examined in the next 
section. 


SIMULATION TEST RESULTS ON A STOL APPROACH TRAJECTORY 


Tests have been carried out using a simulation of the system shown in 
figure 1, including the AWJSRA plant and corresponding acceleration controller 
as well as the Command Generator. Details of the plant and its acceleration 
controller are found in references 2 and 11 but can be omitted from the 
present discussion since the structure of the Command Generator is independent 
of these details. The aircraft dynamics are simulated in a simplified form 
as described in reference 12; that is, the dynamics and the control logic are 
partitioned into translational and rotational degrees of freedom and the 
rotational degrees are replaced with a simple second-order model of attitude 
response to attitude commands while the translational degrees of freedom are 
accurately represented. In addition, the estimation subsystem of figure 1 is 
not simulated; the input signal to the control logic, (fg, Vg, ^g) omits 
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estimation errors due to random measurement noise but simulates the significant 
discontinuities in estimated position that can occur during an approach as a 
result of switching from sensors with large unknown measurement biases (e.g., 
TACAN-DME, barometric altimeter) to more accurate sensors (e.g., microwave 
landing system) . This simulation is well suited to testing the operation of 
the CG and the trajectory control subsystem in response to various reference 
trajectory commands, and to the Trimmap model errors, sensor switching events, 
and external disturbances expected in operational use. 

The reference approach path used for the test is shown in figure 20; it 
is constructed as a sequence of simple legs, and the corresponding parameters 
for input to the CG are included in figure 20. This path includes some 
anticipated features of possible STOL arrival trajectories; compared to CTOL 
operations higher altitude is maintained until close to the runway and then 
descents are steep and confined to limited airspace and with speed reduction 
to the lower STOL landing speeds during descending legs. The test includes 
leg switches requiring multiple axis maneuvering, sensor switch events requir- 
ing unanticipated initialization of the CG from off-reference conditions, and 
both steady and random disturbance inputs (e.g., Trimmap model errors, wind 
turbulence) requiring adaptation of the transition maneuver commands to track- 
ing errors. 

Time histories from the simulation results are shown in figures 21 
and 22. Figure 21 covers the first half of the approach and illustrates the 
capability of the Command Generator to generate acceptable four-dimensional 
maneuver commands, both to execute the known reference approach path and in 
response to unanticipated initializations of the CG to relieve the feedback 
of saturation level tracking errors. Figure 22 covers the second half of the 
path and illustrates the behavior of the CG and the system in the presence 
of both random and steady acceleration disturbances. The variables in both 
figures are grouped into those describing the TD and CG acceleration command 
(e^, e^, aCp in figs. 21(a) and 22(a)), command and response histories 

of velocity coordinates ((V, in figs. 21(b) and 22(b)), the tracking 

errors and corrective feedback acceleration (6rp , 6vp, Aap in figs. 21(c) and 
22(c)), and the applied specific force commands and response ((A^, Afj, «>v) 
in figs. 21(d) and 22(d)). The variables (A^, A^, applied 

specific force supplied by the engine and aerodynamic forces of the aircraft 
in order to fly the commanded trajectory and support the aircraft weight; 
that is , 


Ic " — " (86) 

In cylindrical coordinates and in g units, the specific force along and 
normal to the path (A^ and A]^, respectively) and the orientation of the normal 
force from the vertical plane ($y) , are given in terms of the trajectory and 
path axis acceleration components by 
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here, V_2» P.3} unit vectors defining the path axes reference 

frame. These force coordinates are convenient measures of the implied aircraft 
control activity; the controls are fixed or slowly varying if the applied 
force is constant and otherwise their commanded rates are proportional to the 
commanded specific force rates. For CTOL aircraft, (A^, Aj^, are con- 

trolled nearly independently with throttle, elevator, and ailerons, respec- 
tively. For STOL aircraft, such as the AWJSRA, which may have powered-lift, 
variable-configuration and control redundancy, there can be significant cross- 
coupling of the specific force components among the control variables. These 
details need not be considered in the present discussion since no model of the 
plant dynamics was included in the CG design (G = I) and (A^, A^, ^v) 
therefore be taken as the measure of the control activity required by the 
trajectory commands. It is noted that static equilibrium flight conditions 
(a. = o) and steady turns correspond to fixed values of the applied force 
components. Further, each leg of a reference trajectory in this work corre- 
sponds to a fixed or very nearly fixed value of (A^, Aj^, $v) » complete 

approach path can therefore be viewed as a sequence of such fixed operating 
points with the transition dynamics supplying the transient specific force 
commands required to pass from one operating point to the next. 


Basic Maneuvering Behavior 

The initial 80 sec of the test illustrate the multiaxis maneuvering 
capability of the system in the absence of external disturbances and without 
tracking errors at the time of transition initialization. The reference path 
switches from level flight to a descending turn at t = 11.6 and exits from 
this leg into a level, straight, decelerating leg at t = 45 . 

The TD functions by adding its transition states (e^, e^, e^) to the 
reference states; the initialization of this transition is timed independently 
on each path axis to optimize the maneuver kinematics, and then the three-axis 
maneuver is the superposition of these single-axis transitions. This function- 
ing is illustrated in figure 21(a); the transition is initiated at slightly 
different times for each axis in advance of the reference leg-junction times, 
t^^} , t(^) and the components of the acceleration command, acp, are seen to 
superpose the turn entry and exit on the lateral axis with the pitchover and 
pull-up on the normal axis and a small deceleration step on the longitudinal 
axis. Further, the desired characteristics regarding minimal overshoot of the 
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reference conditions are seen on all three axes in the behavior of , e 2 > 
ac^^), ac(2) in figure 21(a) and (vc, Yc> '^c) figure 21(b). Thus, the CG 
output command retains the desired kinematic features under superposition. 

Tracking performance is given by the tracking errors and corrective 
feedback acceleration shown in figure 21(c). Here, these signals measure 
the success of the CG design structure and parameter values in providing 
suitable feed-forward commands for the given plant-acceleration controller 
system. Ideally, it is desired that these errors be zero in the absence of 
any other disturbance than the CG command since nonzero values imply that 
constraints enforced on the CG command are not necessarily satisfied by the 
aircraft output and that the trajectory regulator is required to use some of 
its margin to execute the reference path. The errors seen in figure 21(c) are 
nonzero temporarily during transitions and are due to control lags excited by 
nonzero acceleration command rates. These rates and, hence, the errors are 
limited through the TD parameter, t. These errors can also be nulled approxi- 
mately if a sufficiently accurate model of the plant dynamics is included in 
the TD, as discussed in the section on maneuver tracking properties. The 
present CG design omits a plant model, but the observed tracking errors and 
feedback acceleration are acceptably small in magnitude and occur only 
transiently. 

Control activity is measured by variation in the applied specific force 
components (A^, A^, defined by equations (85) and shown in figure 21(d) 

for the CG command (^) , the total acceleration command to the Trimmap (f ^|-) , 
and the aircraft response (^) . During a transition, aircraft controls differ 
from their steady-state values for the leg in proportion to the difference 
of the specific force from its steady-state value so that it is desired that 
the transitional specific force commands reach steady-state values without 
unnecessary overshoot, reversal, or activity. However, the command time 
histories of the specific force components can be more complex than those 
of the path axis acceleration components and flight-path angle from which 
they are composed, especially for multiple-axis maneuvering. Referring to 

in figure 21(d), the components ^Vc> ^uc show little overshoot of steady- 
state values and this parallels, respectively, the behavior of ac^^) , 
from the CG. However, A^^,, which is a function of both lateral ana normal 
accelerations, drops transiently at first for the pitchover and then rises to 
its steady-state value for the turn without overshoot, and performs the mirror 
image maneuver at the turn exit. Thus, the required control activity is 
potentially very complex in the case of multiple axis maneuvers but the 
present design economizes this activity as a result of the transition 
initialization algorithm which economizes the acceleration command activity 
from which ^ is composed. The total specific force command, fct> shows 
modestly degraded behavior in figure 21(d) compared to the CG command as a 
result of tracking errors, with higher amplitude excursions from steady state 
and some overshoot, but the specific force response of the aircraft retains 
the desired features. The magnitude of the specific forces required for the 
maneuvers in figure 21(d), remains within ±0.2 g of the static equilibrium 
level flight values, and this range is typical of the modest maneuvering 
domain of passenger operations generally. 
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Response to Sensor Switching Events 

During the time period, [80, 120], the use of the CG to relax large tra- 
jectory errors, which would otherwise saturate the feedback control, and the 
ability of the TD logic to adapt its transition maneuvers to off-reference 
initial conditions is tested. The test consists of a reference trajectory 
turn entry (at t = 96) surrounded by sensor switch events; a switch in 
X— y sensor at 85 sec shifts the estimated position laterally 150 m and at 
105 sec an altitude sensor switch shifts the estimated altitude 20 m. At 
these sensor switches, the CG is signaled automatically to initiate a transi- 
tion back to the reference trajectory. 

The following sequence of events occurs. First, at the x-y sensor switch 
the TD initiates an appropriate transition back to the reference trajectory, 
which in this example is a typical sidestep maneuver (fig. 21(a)). Tracking 
errors develop (fig. 21(c)), principally in the longitudinal direction. 
Secondly, the turn entry of the reference trajectory is initiated in mid 
side step (at 94 sec) from states that are considerably off the reference 
trajectory, both in lateral states, due to the uncompleted sidestep transi- 
tion, and in longitudinal states, due to the tracking errors. The TD adapts 
both the initialization timing and its initial conditions in all three dimen- 
sions to these off-reference conditions (fig. 21(a)); further, the TR is 
relieved by the CG initialization as seen by the discontinuous nulling of all 
error and feedback signals (fig. 21(c)) at this and all transition initializa- 
tions from off-reference states. During the turn entry, longitudinal tracking 
errors again increase. Thirdly, at the altitude sensor switch (t = 105) the 
TD is again initialized from off-reference states and an altitude step 
maneuver is superposed with the turn entry and with the longitudinal tracking 
error relaxation. 

Tracking errors during this period are principally longitudinal and due 
to kinematic cross-coupling neglected in the TD and aroused by lateral axis 
transitions as noted in the previous section. The tracking errors also 
include the effect of unmodeled plant dynamics; that is, transient nonzero 
tracking errors develop in response to the CG command and result in use of 
the TR to execute the commanded trajectory. These effects of this design 
approximation occur only transiently and are of acceptable magnitude in 
the present application. 

The largest transition acceleration commands during this period are in 
the lateral direction. The lateral acceleration and corresponding roll 
commands (figs. 21(a) and (d)) are seen to overshoot their steady-state values 
for the turn; this overshoot does not occur for the reference initial condi- 
tions but cannot be prevented for the unfavorable initial conditions of this 
test, even by the adaptive initialization timing logic. However, the 
acceleration overshoot is limited by the TD excursion limit for the lateral 
axis (0.2 g) and the roll angle command is correspondingly below its 30° limit 
for transient use. 

The longitudinal specific force variations (fig. 21(d)) are dominated by 
off-reference events during this period of the test, first by the longitudinal 
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tracking errors and later by the Y“Variation of the altitude step. The normal 
specific force reflects the combined lateral and longitudinal axis activity; 
the sidestep, turn entry, and the superposition of the altitude step on the 
steady turn can all be detected in its variations. These specific forces 
are complex compared to those of the first test period when off-reference 
events were absent and imply correspondingly complex use of the aircraft 
controls, but the acceleration commands remain within acceptable limits. In 
this test, (corresponding to the CG commands) is piecewise continuous with 

discontinuities at the transition initializations due to off-reference condi- 
tions. The total command, which defines the control commands, is every- 

where continuous as a result of the TD initialization rules. 


Wind Turbulence Effects 


The simulation test includes a period of heavy wind turbulence, 

(135 < t < 195), during which the reference path switches from a decelerating 
level turn to a steep glide slope (at t = 161) and then to a steeper 
decelerating turn (at t = 194) . 

Low altitude wind turbulence is a random vector, assumed generated 

independently along three axes from processes that are stationary for the 
flight duration and Gaussian with variances that depend on altitude and 
direction (ref. 13), but are of the same order of magnitude in all directions. 
Wind turbulence is converted instantly into aerodynamic force variations that 
provide the random acceleration disturbances of the system denoted a^ in 
the system diagram, figure 9. The aerodynamic lift, drag, and side forces 
and their dependence on the air velocity vector can be given as 

F = ^ pva^[-CT)(a)i + CY(g)j - C (a)k ] 

—a 2 ^ ^ —s L — s 


where (i^, jg, kg) are the usual aircraft stability axes, p is the air 
density, (Cq, Cy, Cl) are the drag, side force and lift coefficients, and 
(va, a, 3) are the body axes spherical coordinates of the aircraft velocity 
with respect to the air mass (ref. 13). Denoting the body axis frame as 
(ij 3 , Jl, k^) then (va, a, 3) are defined as 


(va, a, 3) 



and the variation in aerodynamic force due to wind turbulence is 







) 


The first two terms are aligned nearly along since aircraft usually have 

large values of Cl/Cj)* of the order of 10, while the third term is along 
thus, the plant maps the wind turbulence field into acceleration disturbance's 
that are principally normal to the path, as is illustrated by the disturbance 
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time history, (a-ac) , in the present test (fig. 22(c)). The test case dis- 
turbances have large bandwidth and peaks of 0.1 g in the normal direction. 
These result in modest tracking errors and feedback commands, Aa, having 
similar orientation but the lower bandwidth of the acceleration controller 
subsystem. 

The transitions generated by the CG (fig. 22(a)) are little affected by 
the turbulence, even though their initial timing and states are adjusted to 
relieve the TR, because on the longitudinal and lateral axes only small 
disturbances and tracking errors occur while on the normal axis the tracking 
errors are small compared to the discontinuities of the reference trajectory 
pitchovers. The resulting time histories of the trajectory variables 
(ac , vc, ^^) during this period show the desired maneuver command behavior 

with little or no overshoot in lateral acceleration or direction during 
turn exit or entry, or in normal acceleration or flight-path angle during 
pitchovers, or in vc during speed rate changes. It is also apparent in 
figure 22(c) that initializations of the CG relieve the tracking errors due 
to turbulence only momentarily, because this disturbance field is persistent 
in time and the corresponding tracking error statistics are promptly reestab- 
lished. 

The differing character of the specific forces, , ^^9 jfj in the 
presence of turbulence is seen in figure 22(d), especially in the normal 
acceleration component. The CG command, reflects the low frequencies 

of the TD and the constant steady-state values of the reference trajectory 
legs, and is unaffected by turbulence except for discontinuities at TD 
initializations. The total command, superposes the TR feedback on 

and its excursions parallel in amplitude and frequency. Finally, the 

actual force f differs from combined effects of plant dynamics 

and random acceleration disturbance. 


Acceleration Disturbance Effects 

The final test illustrates the effects of acceleration disturbances that 
are constant or slowly varying compared to the characteristic times of the TD, 
and were denoted as a^ in figure 9. These disturbances result from any 
uncompensated mismatch between the actual aircraft forces and the algebraic 
model of aircraft forces used in the Trimmap element of the control system. 
These disturbances can result from large errors in estimating the various 
independent parameters of the force model, such as aerodynamic force coeffi- 
cients, engine thrust model, mean wind, and aircraft weight. 

The effects of such errors is tested in the simulation results by the 
addition of constant drag and lift disturbances of 0.05 g beginning in mid- 
helical turn (at t = 250) . These are significant errors compared to the 
longitudinal acceleration capabilities of the passenger aircraft (usually 
under ±0.2 g) and compared to the desired accuracy for the force model. 

The resulting steady-state behavior of the acceleration controller subsystem 
is seen in figure 22 during the period prior to the glide slope capture; 
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the disturbance is exactly cancelled by a steady feedback acceleration command 
(fig. 22(c)): 


Aa = - ad 


but this equilibrium requires a corresponding position standoff error on the 
longitudinal and normal axes in the amounts: 


6r 


(i) 




i = 1,3 


where are the position error feedback gains. An apparent lateral 

velocity standoff also occurs in the amount, v6r(i)/R^, which reflects the 
steady-state position of the aircraft on the commanded turn but ahead of the 
commanded position by the longitudinal standoff distance, and which is 
nonzero only during turns. 


The above errors are present at the time of glide slope capture (t = 288) , 
so that in the transition commands (fig. 22(a)) the basic turn exit and speed 
rate correction maneuvers required by the reference trajectory are modified to 
relax the longitudinal and normal position and acceleration feedback offsets 
of the TR. However, the acceleration controller feedback promptly restores 
the steady-state offsets, as seen in figure 22(c), by means of a transient 
which roughly cancels that portion of the feed-forward command intended to 
relieve the TR. The possibility of mutual cancellation or reinforcement 
by the two transients is apparent. This process of relieving and restoring the 
TR offsets due to model errors is repeated at each initialization of the TD. 


Further details of the disturbance effects are seen in the specific forces 
(fig. 22(d)). The acceleration commands of the CG (A^^, Ajj^) exhibit 0.05 g 
discontinuities at each initialization of the TD, followed by activity 
generated mostly to relax the TR offsets. The total command, which governs 
the actual control activity, is better behaved, and shows both continuity 
and less activity as a result of the mutual cancellation of feedback and TD 
transients. The plant output also maintains these favorable properties and 
settles in steady state to the CG acceleration command as desired. 


In these tests the steady disturbance was compensated by the feedback 
control; the resulting offsets of 6r and Aa interact with the adaptive 
initialization of the transition dynamics to produce ineffective transition 
commands from the CG on the longitudinal and normal axes as noted above. 

These can be removed, along with the offsets in by providing compen- 

sation for ^ in the Trimmap element of the control system, as noted in 
reference 2. Any residual uncompensated disturbance would, however, result in 
the interaction between CG and the tracking system at leg initialization times 
noted above and, if of a significant size, further disturbance estimation and 
compensation in the Command Generator would be useful. 
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DISCUSSION AND CONCLUSIONS 


This report has described the development of a trajectory Command Genera- 
tor for a class of advanced digital flight-control systems (ref. 1; fig. 1) 
with capabilities of executing complex four-dimensional terminal area flight 
paths as directed by an air traffic control and guidance system (ref. 3) or 
selected by a pilot. The design methodology was presented, including the 
objectives, structure, and a detailed synthesis procedure for the CG. The 
CG described here was further examined in an autopilot of the type defined in 
reference 1 applied to the AWJSRA aircraft and subjected to simulation testing. 
Factors studied include the nature of the basic operational repertoire of 
maneuvers and operational constraints, and the system behavior in response to 
four-dimensional maneuver commands and as affected by various CG design 
approximations and disturbances found in the operational environment. These 
simulation results indicate that a satisfactory nonlinear system was developed 
which satisfies the basic design objectives while maintaining a practicable 
degree of simplicity. Further, since the CG structure and many of the opera- 
tional constraints are independent of the aircraft, the simulation results are 
representative of STOL and CTOL aircraft and of passenger operations generally. 

The CG design problem in the context of the system of figure 1 is to 
generate suitable feed-forward trajectory commands to a trajectory tracking 
subsystem in order to execute an input reference trajectory. The class of 
input reference trajectories assumed for the problem is typical of passenger 
operations, including those of advanced terminal area ATC systems. These are 
sequences of four-dimensional segments that can be straight lines with constant 
inertial acceleration or circular arcs with constant radius each of which 
satisfies operational constraints on steady-state flight. The task of the CG 
is to provide maneuvers that transition the aircraft from one steady flight 
condition to the next and satisfy the design objectives for such maneuvers. 

The input sequence also possesses certain approximate properties in practice 
which influence the design; that is, the duration of each segment is usually 
longer than the time required to transfer from the previous segment and the 
transfer between successive segments can usually be made optimally in the 
vicinity of the segment junction times. The design criteria relating to the 
output are that the aircraft trajectory satisfies operational constraints 
on steady-state and transient kinematics, {gi(^, A) i 

and that transition maneuvering be optimized to reduce or minimize accelera- 
tion activity and overshoot of the steady-state flight conditions; these are 
considered necessary conditions for operational acceptability of the system. 
While these operational constraints are to be satisfied by the actual aircraft 
trajectory it suffices to enforce them on the feed-forward trajectory commands 
of the CG provided acceptable tracking errors are obtained from the given 
aircraft and trajectory tracking system over the desired operational domain 
of the system. 

The CG function is carried out in the proposed design using two principal 
elements: the transition dynamics (TD) and its initialization algorithm. 

The trajectory discontinuities between successive legs are viewed as initial 
errors from the new segment which the TD is to relax to zero in accordance 
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with the desired operational constraints. A TD was synthesized to generate 
these transitions independently on each of the aircraft path axes and a 
dual-mode control law was synthesized for each axis which is globally asymp- 
totically stable with (1) selectable damping and characteristic time of the 
linear mode transients and a maximal region of initial conditions using linear 
control; (2) a nonlinear mode with selectable parameters to limit velocity 
and acceleration excursions about the reference flight condition; and (3) a 
selectable limit on the maneuver acceleration rates in order to enforce band- 
width restrictions on the maneuver commands. System behavior with this TD 
was examined in the third section of this paper. Initial conditions for the 
transition were derived to null initial tracking errors and maintain contin- 
uity of the acceleration command to the plant. Further, it was found that 
perfect tracking was obtained if disturbances were constant and exact models 
of the plant dynamics and disturbances were incorporated in the TD. The TD 
structure and control law synthesis were extended to include these models, 
G(s), a"d. The potential value of minimizing tracking errors is that (1) the 
use of the trajectory regulator to execute the feed-forward commands is 
minimized so that its control margin is reserved for the regulation of random 
unpredictable disturbances; (2) the kinematic constraints are more accurately 
imposed on the plant output; and (3) the bandwidth of the TD commands for 
which acceptable system behavior can be obtained is increased. The applica- 
tions work of this report omitted such models (equivalently, the design 
assumes G(s) =1, ad = 0) but satisfactory results can be expected with this 
design provided (1) the plant dynamics are of sufficiently higher frequency 
than the transition dynamics and (2) the steady disturbances are sufficiently 
small. The transition initialization switching logic (established in the 
fourth section) provides optimization of the transition maneuver. The 
proposed optimization criteria were to minimize initial jerk for transitions 
initiated in the linear mode, and to prevent saturated control reversal and 
corresponding peak jerk in the nonlinear mode. In addition to reducing the 
local acceleration activity specified by the criteria, transitions are 
obtained with approximately minimum duration and overshoot of the reference 
flight conditions as was shown in extensive examples. 

The scope of aircraft operations to which this design is applicable is 
bounded in several ways. First, the problem formulation assumes a class of 
input reference trajectories typical of passenger operations. The properties 
of this class are reflected in the CG design both in the concept of transition 
dynamics and in the use of segment initialization time criteria. However, 
the flight-control structure of reference 1 is applicable to more general 
aircraft operations and a larger class of inputs is generated by such opera- 
tions than can be accommodated by the present CG design. Secondly, the gener- 
ation of independent transitions in path axes neglects Coriolis accelerations 
which cross-couple the excursion states during transitions with significant 
path axis angular velocity. This approximation is not expected to remain 
adequate for operations outside the modest maneuver domain and angular 
velocities of the passenger operations postulated in this study. Thirdly, 
the path axes frame is undefined at special points in the domain of helicopter 
and VTOL operations (hover and vertical flight) and an alternative choice of 
axes in which to generate transitions is appropriate for these classes of 
aircraft. 
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An algorithm for simulation test and study was presented; it has a 
simple structure and utilizes equations directly from the text. The class of 
input reference trajectories and its properties were examined in detail and 
equations for an appropriate continuous- time reference trajectory generator 
were given. The period over which the TD could be initiated by the switching 
algorithm was limited to about 10 sec in advance of the reference leg— junction 
times based on properties of the input set. Further, it was noted that the 
desired values of parameters constraining then transition excursions were 
functions of the reference flight condition and the algorithm provides for 
their evaluation at the beginning of each leg. Use of a detailed model of 
plant dynamics, G(s), in the TD was not included in the present simulation 
study in view of lower bandwidth of the trajectory commands compared to plant 
dynamic s . 

Parameter values for the TD were established for the present study on 
the basis of operational constraints on trajectory kinematics. In general, 
these constraints are an essential part of the design problem associated with 
the operations for which the design is to be acceptable, especially during 
approach and landing when the aircraft operates near the boundaries of its 
performance capabilities. These constraints require a CG structure that 
recognizes the nonlinearity of the problem and provides for their systematic 
enforcement. In the CG design problem, some constraints are applicable to 
the class of operations (passenger comfort constraints) and some are derived 
for the specific aircraft from safety margin constraints and aircraft perfor- 
mance limits. The appropriate saturation parameters for the TD can then be 
derived from these constraints as functions of the reference flight condition. 
The basic repertoire of single-axis maneuvers was also examined to illustrate 
the significant kinematic variables in passenger operations and to demonstrate 
the general effectiveness of the TD saturation parameters in limiting excur- 
sions and of the TD initialization switching logic in optimizing these 
variables. The effects of Coriolis accelerations neglected in synthesizing 
the transitions was also examined; these result in kinematic incompatibility 
of the trajectory commands and are significant only transiently as longitudinal 
tracking errors during lateral axis transitions. These errors are proportional 
to initial values of the lateral TD states, but can be limited when large side- 
step distances or large direction steps are to be executed by elaborating the 
input reference trajectory. 

Finally, behavior of the CG and the flight-control system was tested on 
a STOL approach test trajectory and in response to the principal disturbances 
present in operations. These tests indicate a number of conclusions. 

First, for the postulated class of input trajectories the CG can 
synthesize satisfactory four-dimensional, multiple-axis transitions from the 
single-axis transitions. The optimal timing and excursion saturation limits 
imposed independently on each axis remain effective in the multiple-axis 
command in minimizing overshoot and limiting excursions and rates of the 
kinematic variables of interest and of the corresponding aircraft control 
activity given by the specific force commands. 
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Secondly, the transition dynamics adapt successfully to off-reference 
initial conditions of significant magnitude due either to unanticipated 
initializations as a result, for example, of sensor switch events, or of 
uncompleted earlier transitions. In the latter case, the optimal timing logic 
also adapts to the off-reference condition. The bounded set of initial TD 
conditions over which the system functions acceptably and for which these two 
conclusions are operat?lonally valid cannot be established in the present 
simulation. However, the simulation test path from which these conclusions 
are drawn contains a rigorous set of example maneuvers of greater complexity 
than is expected operationally. 

Third, unmodeled or inaccurately modeled plant dynamics in the TD intro- 
duce tracking errors during all transitions and involve the regulator in 
executing the reference trajectory. Similar effects occur from unmodeled 
Coriolis accelerations as noted earlier. The effect on performance of these 
simplifying approximations were found acceptable in the simulation tests but 
of sufficient magnitude that further development of the theory and its 
application is of interest in both cases. 

Fourth, the CG adapts the initial conditions for transitions to any 
steady-state or time-varying tracking errors of the acceleration controller 
subsystem. If these tracking errors are due to persistent steady disturbances 
(e.g., uncompensated Trimmap model errors) or to random disturbances with 
persistent statistics (e.g., wind turbulence) then the trajectory tracking 
subsystem necessarily restores the same steady-state error or its statistics 
that prevailed before the transition. Thus, ineffective transient activity 
is added to the transition commands by this adaptation. The effect is 
proportional to the magnitude of the tracking errors and is acceptable for 
sufficiently small disturbances. The effects of wind turbulence were found to 
be small while steady disturbances can be compensated within the trajectory 
tracking system. Persistent large disturbances which saturate such compensa- 
tion cannot be effectively relieved by the CG and some other remedy, such as 
revising the reference flight condition, is indicated. However, the ineffec- 
tive transient command activity due to uncompensated steady disturbances 
can be minimized by inclusion of disturbance estimation and compensation in 
the Command Generator. This problem is also of interest for further study. 

Continued simulation study of the proposed CG algorithm is planned using 
increasingly realistic simulations of the flight-control system and aircraft, 
and including the estimation and attitude control subsystems. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, California 94035, May 2, 1978 
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APPENDIX A 


DERIVATION OF DUAL-MODE CONTROL LOGIC 


General Solution for the Region of Linear Contrbl 
Consider the system 


X = Ax + Bu (Al) 

where x is an n-vector, u an m-vector, and A; B are constant matrices 
of appropriate dimensions. Then the control law u that minimizes the 
performance measure 

oo 

Q (x^Qx + u"^Ru)dt (A2) 

is given by reference 14 

u = -k'^x (A3) 

where 

K = PBR'^ (A4) 

and P is the positive definite, symmetric solution of the Riccati equation 

PA + a'^P - PBR~1b^P + Q = 0 (A5) 

and it is assumed that R > 0 and (A, Q%) is completely observable. For this 
control law, the optimal cost function is 

V = x^Px (A6) 

which is a Lyapunov function for (Al) and (A3) with the time derivative 

V = -x'^Qx - x^PBR“^B^Px (A7) 

Suppose now that the control law u(x) is subject to saturation in the 
following way: 



-Urn,- ^ u^(x) < 

X 

(A8) 

where the subscript (i) refers to the ith component of a vector, 
therefore permit the control to be given by 

We can 


u . (x) = -k.*^x 

X X 

(A9) 

for those 

X satisfying ^ 

< -k^ X < u^^ 

(AlO) 
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Here, refers to the ith column of K. Let us now further restrict u by 

finding a region C so that if 

x(to)^^ (All) 


then 

x(t) E 0, and u = -K^x , t > t^ 

that ±Sy Q is a region of asymptotic stability as well as a region in which 
the control is entirely linear; Q is an invariant region of asymptotioally 
stable linear control. 

A simple and convenient way of finding one such region, is offered 

by the function V(x) in (A6) . Let Xj_ be the two lines defined by 





(A12) 


and then compute 


= min I min [V(x)]| (A13) 

l<i<m 

which is found to be (cf. the case m = 1 in refs. 15 and 16) 

l<i<m ‘ ^ ^ 


where 

Then: 

Qq = (A15) 

Now, a dual-mode control scheme may be constructed based on inside 

Qq use the law (A3) and outside use a nonlinear, saturating control law. 

Such a controller will have the property that (ignoring disturbances) , once 
the linear region is entered, the control remains linear. Kuznetsov (ref. 4) 
was apparently first to use such a method (for m = 1) to construct dual-mode 
controllers. He essentially used a surface like that in (A15) to define the 
boundary between control regions and explored various possible constructions 
of the nonlinear controller. The next section presents the details of a 
particular dual-mode controller for the transition dynamics of this paper. 

An invariant linear region larger than Qq is derived and the control law 
for the nonlinear mode is given. 


u = mini 

o . ( 
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Dual-Mode Control Law for m = 1, n = 2 


The invariant linear region obtained in (A13) is generally only a 

subset of the largest such region. Since we put a premium on linear control, 
it would be desirable to be able to enlarge this region. For the case of 
interest here (m = 1), Qq can be extended by various methods, including direct 
numerical calculation of the largest such region as was done, for example, 
in reference 6. We adapt the work of Willems (ref. 5) to this purpose because 
of the convenience of the required calculations.^ 


Consider then the case m - 1. Using Willems’ results, it can be shown 
that the region Q: 

Q = |x:-Uj^ < k*^x < u^ ; x'^Px ^ ^ (A16) 


is an invariant region of asymptotic stability for the linearly controlled 
system, where 


V’ 

M 


Vm 


A 


k^AP ^A\[k'^P ^kk'^AP - (k\p ^k)^] 

A - Bk'^ 


V’ 




M 






(A17) 


(A18) 


and the equality in (A18) holds only in special cases. 

A simple control law for x^f2, which satisfies the saturation constraint 
(A8) , is the usual time-optimal control (i.e., the control that would minimize 
the time to the origin in the absence of the linear control mode in the 
region . This control provides the switching logic shown in figure 6(a) 
excluding the modifications for velocity limiting. This control also has the 
virtue of approximately minimizing the time spent outside The minimum 

time control to the region Q may also be found; it is generally different 
from the minimum time control to the origin, although the increased complexity 
does not justify its use in our application. 


^Our basic method of constructing a dual-mode control is not essentially 
tied to a quadratic optimal law; any as 3 miptotically stable linear control law 
could be used, with any one of an available set of quadratic Lyapunov func- 
tions used to define the region. The approach here is felt to be somewhat 
more motivated, however, since it gives rise to a region that has significance 
in terms of a measure of control performance. 
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Specific Control Law for Equation (22) Without Velocity Limiting 


Considerable simplification is achieved in the implementation of the 
dual-mode control law when applied to the system described by equations (22a) , 
that is, the system 



(22a) 


Without loss of generality, we take R = 1- It can also be shown that 
does not affect the locations of the characteristic roots 

controlled system in this case; hence for simplicity we take ^ Then, 

expressing all results in terms of the easily understood parameters 
instead of the less intuitively appealing parameters Qii, Q22» obtain 


^^2 = 


-2 


(A19) 


where the characteristic roots are now 


^ 1,^2 = ("^ - - l |) 


and 


^12 


-2 


P = 2t ^ 
^22 ^^L 


(A20) 


P = P P 
^11 "^12 22 


and, after some calculation, 

Vm = T ?2(3254 _ 24^2 + 6) 

o ^ 


(A21) 


These results all require that2 




> 


1 

2 


2Although it is unlikely to be desired, characteristic roots with 4 < 1/2 
could be obtained by procedures different from the optimal quadratic synthesis 
above . 
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by 


The minimum- time control for equation (22a) and for x ^ is now given 


where 


u 


a 


f 

{ 


■u 

m 

Q 

V 

O 

or (a 

= 0, 

X > 

0) 


Q 

A 

O 

or (a 

= 0, 

X < 
2 

0) 


U X 


o 

A 

m 1 

2 2 

z 


1 2 
2 ^2 

x^ < 0 
2 


(A22) 


(A23) 
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APPENDIX B 


VELOCITY-LIMITING MODIFICATION OF CONTROL ALGORITHM 


Appendix A described a dual-mode control algorithm that maintains 
prescribed saturation limits on the scalar control u in the form 


u 



for 


V^(A22) for 


X 

X ^ ^ 


(Bl) 


where Q is defined by equation (A16) . Here the law is modified (for n = 2) 
to maintain, in a certain sense, bounds on the velocity 


-X 


2 m 


< X 


2M 


(B2) 


The modifications are developed next by considering different types of initial 
conditions. 

In the case that x(0) is in ^ and x^(0) satisfies equation (B2) , a 
review of phase plane trajectories and the basic control algorithm equa- 
tion (A22) shows that a trajectory either enters or else encounters a 

limit of equation (B2) , say, X 2 = (see fig. 23). If equation (Bl) is 

modified so that u = 0 at this point, X 2 will remain constant at the 
limit value, and this control should be held until the first opportunity to 
decrease X 2 by reinstating the basic algorithm (eq. (Bl)). This can occur 
at either of two events; first, if the switching curve a = 0 is encountered 
before the boundary of Q (trajectory i^ in fig. 23), then the control 
simply switches to u = -Ujj^ at the switching curve in accordance with 
equation (A22) . Secondly, if the boundary of 0. is encountered before the 
curve a = 0 , then the control remains null but switches to the linear control 
law of the region 0. as given by equation (Bl) when k*^x = 0 (see fig. 24). 

We need only demonstrate that the modified linear region as shown in 

figure 24 (region U with the "bite" taken out of it) as an 

invariant one, that is, that trajectories cannot cross the new boundary-line 
segments 0 x *“02 and 02-03 from interior points of This is enforced for 

0 ^- 02 , a locus of constant X 2 , by the nulled control of the modified law 
along this segment. For segment O 2 -O 3 , observe that x^ is strictly 
increasing: 

x^ < 0, x^ = x^ > 0 


and the trajectory is locally parallel to the x^ axis: 


dx. 

dx. 


1 T 

-k X 

Xo 


= 0 
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so that a trajectory necessarily evolves into the interior of from any 

point on 02-03. Since the linear control becomes negative at interior points 
of adjacent to O2-O3, the velocity subsequently decreases. To conclude 

this case, note that similar arguments apply to maintaining the limit ~'^2m 
in the fourth quadrant of the phase plane. 

For initial conditions with x( 0 ) in Q but X2(0) does not satisfy 
equation (B2) , there are two subcases to consider; first, if 

x^CO) < and a < 0 

or 

x^(0) > x^j^ and a > 0 

the control law (eq. (Bl)) need not be modified since it already drives the 
velocity toward the permitted range (eq. (B 2 )) as rapidly as possible. The 
second subcase occurs with 



" -*2m 

and 

Q 

V 

0 

=■ “‘ZM 

and 

Q 

A 

0 


(B 4 ) 


Here we can consider modifying the control to reduce x^ at the maximum rate 
as for the first subcase, but the possibility of *'chatter**-like phenomena 
whenever the velocity limit line (where u = 0) is encountered must be 
accounted for. Alternatively, we select u = 0 for the entire region given 
by equations (B 4 ) ; the resulting free trajectory must exit this region in 
finite time. 


Finally, the remaining possible initial conditions, x( 0 ) G Q either 
begin in the ”bite” removed from the linear region, in which case the control 
is again nulled as it is on the region defined by equations (B 4 ) , or it 
begins in where the linear control law of equation (Bl) remains in force. 

In summary, equation (Bl) is simply modified to give null control on the 
following regions: 


for 

or 


u(x) = 0 
T 

X G and { (“k x < 0,x^ < -x^jj^) 
K G and { (a > 0,x^ < x^^) 


or 


or 


(-k^x > OjX^ > 
(a < 0,X2 > 


(B 5 ) 
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Figure 1.- Structure for an automatic flight-control system. 



Figure 2.- Trajectory command generator structure 



















(a) Control law and switching boundaries in the x-plane 
Figure 6.- Transition dynamics. 



Figure 6.- Continued. 
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(c) Xp x^ plane trajectories — optimum linear 
transition. 
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Coiranand generator algorithm 
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(a) Computational flow diagram. 

Figure 14.- Transition initialization algorithm. 
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(b) Summary of equations for transition initialization algorithm. 


Figure 14.- Concluded 
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(a) Computational flow diagram. 

Figure 15.- Transition dynamics algorithm for simulation test. 
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REQUIRED INPUT PARAMETERS FOR TD: 
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(b) Summary of transition dynamics equations 
Figure 15.- Continued. 
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(b) Concluded. 


Figure 15 


Concluded 




Figure 16.- Lateral axis 


single state transitions 
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(b) Effects of transition initialization time 
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(c) Longitudinal axis activity 


Concluded. 
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(b) Effects of transition initializa- 
tion time. 


Figure 18,- Turn entry transition. 
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(a) x-plane control regions and transition 
trajectories* 


(b) Position step transition. 



(c) Velocity step transition. 


Figure 19.- Longitudinal axis transitions. 
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Figure 20 


STOL simulation test approach path 







(a) Command generator variables. 

Figure 21.- Simulation results — first half of test path. 
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(b) Velocity commands and responses 
Figure 21.- Continued. 
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(c) Trajectory tracking variables. 
Figure 21. ~ Continued. 
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(d) Applied specific force commands and responses 


Figure 21.- Concluded. 
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(a) Command generator variables. 

Figure 22.- Simulation results — second half of test path. 
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Figure 22.- Continued. 
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(d) Applied specific force commands and response. 
Figure 22,- Concluded. 
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